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A^ ' We present a new method for the evaluation of the perturbative expansion of the QCD 

V^ , pressure which is valid at all values of the temperature and quark chemical potentials in the 

(^ ' deconfined phase and which we work out up to and including order g^ accuracy. Our calcu- 

xl , lation is manifestly four-dimensional and purely diagrammatic — and thus independent of 

any effective theory descriptions of high temperature or high density QCD. In various limits, 

we recover the known results of dimensional reduction and the HDL and HTL resummation 

schemes, as well as the equation of state of zero-temperature quark matter, thereby verifying 

their respective validity. To demonstrate the overlap of the various regimes, we furthermore 

show how the predictions of dimensional reduction and HDL resummed perturbation theory 

Oh. agree in the regime T ^ ^/gn. At parametrically smaller temperatures T ~ g/i, we find that 

the dimensional reduction result agrees well with those of the nonstatic resummations down 

to the remarkably low value T « 0.2mD, where r?iD is the Debye mass at T = 0. Beyond this, 

we see that only the latter methods connect smoothly to the T = result of Freedman and 

k> , McLerran, to which the leading small-T corrections are given by the so-called non-Fermi- 

^ • liquid terms, first obtained through HDL resummations. Finally, we outline the extension of 

5t I our method to the next order, where it would include terms for the low-temperature entropy 

and specific heats that are unknown at present. 
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I. INTRODUCTION 

The most fundamental thermodynamic quantity in the theory of strong interactions, the QCD 
pressure Pqcu{T, n), can at large values of the temperature T or the quark chemical potentials 
/i be computed in a weak coupling expansion in the gauge coupling constant g, defined in the 
MS renormalization scheme. In the region where T is larger than all other relevant mass scales 
in the problem, the expansion has been extended to include terms of order ^r^log^ [l|-|al) while 
at T = and /x much greater than the critical chemical potential Hc, the pressure is known up 
to and including terms of order g^ 9]. In between these regimes, at < T < gfi, anomalous 
contributions from non- Fermi- liquid behavior have been obtained [li|, which involve fractional 
powers and logarithms of g. The purpose of this paper is to connect all these disjoint computations 
through one expression which gives the pressure at all values of T and /x up to and including terms 
of order g^. 

Ultimately, the reason for the existing results being valid only in the separate domains de- 
scribed above is that different computational methods are practical in different regions of the ^-T 
plane. When T is larger than all other dynamical scales, one can, in the spirit of effective the- 
ories, integrate out the degrees of freedom corresponding to non-zero Matsubara modes (having 
thermal masses of order ttT or higher) to obtain a simpler three-dimensional effective theory. For- 
mally, the requirement for this is that one must have^ T 3> mn which explicitly excludes the 
region of T < gfi. At these parametrically lower temperatures, perturbation theory requires a 
different reorganization which can be most efficiently performed via the hard dense loop (HDL) 
approximation. This approach corresponds to a different effective field theory which is intrinsically 
four-dimensional and whose non-local effective action, known only to one-loop order, was found 
by Braaten, Pisarski, Taylor and Wong [1^. It involves a resummation of nonstatic self-energies 
which, in fact, was a crucial ingredient already in the classic calculation of Freedman and McLer- 
ran |9| for the 0{g'^) zero-temperature pressure. In practice, even before reaching this limit, at 
T ~ ij.g~^ exjp[—37r'^ / {^/2g)] [ij], one of course encounters the non-perturbative pairing instability 
(color superconductivity), but in a strictly perturbative calculation these effects can never be seen. 

In view of the fragmented status of the various perturbative results on the fi-T plane, there is 
obvious motivation for attempting to develop an independent and uniform method of calculation 
which would have the power of both verifying the validity of all the existing computations and 
providing a smooth interpolation between them. Important steps in this direction have already been 
taken through recent advances in the analytic calculation of sum-integrals at arbitrary temperatures 
and densities |3, llJ] as well as in the numerical evaluation of multi-dimensional integrals involving 
the one-loop gluon polarization tensor [l5|. In fact, these techniques provide all the required 
machinery for the first purely diagrammatic four-dimensional determination of the QCD pressure 
up to order g"^ at arbitrary fj, and T which is what we set out to perform in this work. We take as 
our starting point the systematic analysis of all relevant classes of Feynman graphs contributing 
to the partition function up to and including this order which implies that our calculation will be 
independent of any effective descriptions of the fundamental theory that are valid only in limited 



Due to the chemical potential ^ appearing in the fermionic propagator with an imaginary unit, this parameter 
contributes to the long-distance behavior of the free fermionic correlation function (i.e. the Fourier transform of 
the free propagator) only through an oscillatory phase factor. Therefore, from the point of view of IR physics it can 
be identified as an irrelevant parameter, whose value only affects the validity of dimensional reduction through the 
Debye mass scale. Note, however, that this reasoning does not apply to imaginary values of the chemical potential, 
in particular when its magnitude becomes comparable to ttT. (In this region one is in any case restricted to values 
/i < nT/Nc due to the loss of periodicity in the pressure as a function of imaginary fi in standard perturbative 
calculations |ll| . This is due to the explicit breaking of the Z{Nc) symmetry, whose effects are otherwise partially 
visible even in the presence of dynamical fermions.) 



regions of the fi-T plane. 

As is usual in a weak-coupling analysis, we regard the coupling g as an in principle arbitrarily 
small parameter, in which we set out to expand the pressure around the value g = 0. In practice, 
the QCD coupling scale is of course about 1 at 100 GeV and at 1 GeV even about 2, so for our 
results to be applicable we must assume to be working at sufficiently high temperatures and/or 
densities. A critical region where the nature of the perturbative expansion of the pressure changes 
both qualitatively and quantitatively turns out to be that of parametrically small T < rrio ~ g^i, 
so we have found it useful to quantify the smallness of T by expressing it in the form T ~ g^fJ-- We 
then study the changes that occur when x is increased from zero to values larger than 1, where the 
pressure can be re-expanded in a series of (fractional) powers and logarithms of g. 

Our main results, displayed in Sec. IIV| include the following. In the region T//i ~ g^, we 
explicitly recover the perturbative expansion of the pressure to order g"^ which was derived in 
Ref. [a] using dimensional reduction (DR). We show that this result is valid up to arbitrarily high 
values of the ratio /x/T, verifying the non-trivial behavior of the different perturbative coefficients 
as functions of the parameter. In deriving these results, we notice that the larger we choose fi/T to 
be, the smaller values of g we have to use in the series expansion of our numerical result in order to 
find agreement with the dimensional reduction one. We interpret this as a demonstration of the fact 
that dimensional reduction amounts to expanding the pressure in powers of mo/r on top of g and is 
thus inapplicable in the region T < g//. For temperatures T ~ ^/g|J', which still fall into the domain 
of validity of dimensional reduction, we demonstrate the overlap of the different existing results by 
deriving the same expansion for the pressure (in half-integer powers of g) from both the perturbative 
results of Ref. |fl| and the HDL resummation scheme of Refs. jlfll.llfil|. At parametrically even smaller 
temperatures, T < gfj,, we then show how our new calculation — the first one entirely independent 
of the original one of Freedman and McLerran [2| — smoothly approaches their famous result for 
the T = pressure, at the same time reproducing the "anomalous specific heat" contributions to the 
pressure obtained previously through HDL resummations. In contrast, the dimensional reduction 
result is observed to agree numerically remarkably well with those of the nonstatic resummations 
down to T ~ 0.2?tid, or 27rT ~ 1.3rn-D, but it diverges logarithmically in the limit T — > 0. We 
explain how this divergence is related to the behavior of the physical, resummed expression for the 
pressure at large but order g^ values of T/{g^). 

The paper is organized as follows. In Section ^ we review the most important methods ap- 
plied previously to the evaluation of the QCD pressure, namely dimensional reduction, HTL/HDL 
resummations and zero-temperature perturbation theory. After this, we present our new scheme 
in Section II I II where we identify the classes of diagrams we need to compute and then write down 
our new expression for the pressure as the sum of an analytic power series in g'^ and a logarithmic 
sum-integral that is to be treated numerically. In Section HVl we then finally display and discuss 
our results. We plot them in various regions of the /x-T plane as characterized by the exponent x 
in T ~ g^'/x and compare our predictions to the existing results of dimensional reduction ;§] and 
HTL/HDL resummations [ifll. Ilfi|. After this, we summarize our findings in Section where we 
explain in a detailed fashion the structure of the perturbative expansion of the QCD pressure on 
the n-T plane, and then draw final conclusions and briefly look into new directions in Section IVll 
Several computational details as well as the explicit analytic forms of many individual contributions 
to our result are finally displayed in the Appendices where, in particular. Appendix ^ contains 
details on the numerical evaluation of the logarithmic sum-integral introduced in Sec. II I II 



A. The notation 

To conclude the introduction, let us fix our notational conventions which to a large extent follow 
those of Ref. [3]- First, throughout the text we will be working in the Feynman gauge, in which 
many of the existing results our calculation relies on have been derived. The chemical potentials 
of the different flavors, for which we apply the notation 

^^2^' ^ = 1/2-^/^' (1-1) 

are assumed to be equal, and thus we suppress all subscripts for them. (The renormalization scale 
will be denoted by the symbol A and should not be confused with the above p..) For the analytic 
result of Eq. ()3.2U|) . the generalization to several independent chemical potentials is trivial and only 
requires the multiplication of the results by 1/Nf ^ r, but for the outcome of the numerical calcu- 
lation, this is less straightforward. All of our methods can naturally be applied to any combination 
of T and the /i/'s, but this will have to be done on a case-by-case basis. 
As usual, we denote the most common group theory constants by 

CaS""^ = /"''^/"^^ = Nd^'^, (1.2) 

CpSij = {T°-T°-)ij = ^^ Sij, (1.3) 

Tpd"^ = TrT'^T'' = ^<5"^ (1.4) 

dA = 5'"' = iV^ - 1, (1.5) 

dp = ^ii = dATp/Cp = Nnf, (1.6) 

where the T^ stand for the canonically normalized generators of SU(A^) and the trace is taken over 
both color and fiavor indices. 

Our notation for the special functions that appear in the e-expansions of sum-integrals at finite 
density follows that of Ref. |a| , which introduced 

C{x,y) = d^C{x,y), (1.7) 

^n,z) ^ (:>{-n,z) + {-ir+\'{-n,z*), (1-8) 

H(z) = ^(z) + ^(z*). (1.9) 

Here, ^ stands for the generalized Riemann zeta function and ^ for the digamma function 

*H . ^^. (1.10) 

Finally, we will throughout the following chapters apply dimensional regularization at d = 4 — 2e 
dimensions and (unless otherwise stated) in Euclidean metric, using as the momentum integration 
measure 

f ^ f d'^'^p ^ 2e f d^-'p ^ /eU^Y f d'^-^p 
Jp - J (27r)rf-i J (27r)rf-i \ ^tt J J (2^^-^ ' ^''''^ 

where A = (47re~''')^'^A is the scale parameter of the MS scheme. The symbol P/{P} refers to 
bosonic and fermionic momenta, respectively, for which capital letters stand for four-momenta and 



lower-case ones for three-momenta. Using these definitions, bosonic and fermionic 1-loop master 
integrals are defined by 

The bosonic and fermionic distribution functions at temperature T and chemical potential ji are 
denoted by 

nm = ^^^, (1.14) 

"/(^) = \ ( e(^-M)/^ + l + e(^+/^)/^ + l ) ■ ^^-^^^ 

II. PREVIOUS METHODS 

Let us begin by reviewing the main results available for the QCD pressure through various 
previous perturbative calculations in order to summarize the current understanding of the behavior 
of the quantity in various regions of the ^i-T plane. The methods to be covered include dimensional 
reduction which has been used to determine the pressure up to and including order g^lng at finite 
T and /^ |3) la| ) the perturbative T = techniques of Freedman and McLerran [9| , and the HDL 
resummations that have proven to be highly effective at temperatures non-zero but parametrically 
smaller than /x (lO|, [l^] . 

Recently there has been some progress to include effects of nonzero quark masses in the dimen- 
sional reduction approach [lg| as well as to the next-to-leading order pressure at zero temperature 
[]jj. In the following we shall always restrict our attention to the case of massless quarks. 

A. Dimensional reduction 

In finite-temperature QCD, dimensional reduction is based on the observation that at suffi- 
ciently high temperatures one has the parametric scale hierarchy m^^^ <C rrio <C i^T between the 
magnetic mass ~ g^T, the Debye mass ~ gT, and the smallest non-zero Matsubara frequency, 
respectively. One may then integrate out all the non-static degrees of freedom from the original 
theory (corresponding to the last of the above energy scales), leaving behind a three-dimensional 
effective theory for the static gluon fields. At leading order, this integration-out can be performed 
explicitly, but beyond that it becomes more convenient to start from the most generic possible 
three-dimensional Lagrangian that respects the correct symmetries and to determine the values of 
its parameters by matching a set of simple physical quantities to their values in the full theory. 
This leads to a Lagrangian for the effective theory, commonly dubbed electrostatic QCD (EQCD), 
of the (Euclidean) form J>, JO] 

/^EQCD = ^ TV 4 + Tr [( A^o)'] + m|Tr AI 
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o3 



+ f2E^/Tr^o + Ai'U^^o)' + A^'^TY4 + <5£,. (2.1) 

In the above Lagrangian, m^ agrees to leading order with the physical Debye mass mo. 



3 TT"^ 



"^D = g\l ": T2 + Tpf^ , (2.2) 



and the last term 5C^ corresponds to higher-dimensional operators that can be neglected for most 
practical purposes, as they contribute e.g. to the pressure starting only at 0{g'^). The cubic, C-odd 
term [20|, |21i] in Eq. H2.1|) is only present at finite chemical potentials and contributes to the pressure 
at order g^hig, (but turns out to be responsible for the leading order result of the off-diagonal quark 
number susceptibilities at zero chemical potential |2j]). The fields have been rescaled from their 
four-dimensional values by \/T to have the canonical dimensionality at d = 3, and the field strength 
and covariant derivative contain the three-dimensional gauge coupling g^ = g\/T + 0{g^). The 
new theory should be able to describe the physics of the original one at distance scales of order 
\/{gT) and higher. 

Assuming g to be sufficiently small {i.e. T sufficiently high), the above integrating-out procedure 
can be continued one step further by removing also the massive longitudinal gluon Aq from the 
theory, thus producing an effective three-dimensional pure Yang-Mills theory for the Ai fields. This 
theory, magnetostatic QCD (MQCD), is defined by the Lagrangian 

with the gauge coupling constant ^m = 5e + 0{g'^/mE) appearing in the field strength tensor. The 
perturbative expansion of the pressure of the full theory can then be given in the form 

p = Pe+Pm+Pg, (2.4) 

where the first term corresponds to the coefficient of the unit operator of EQCD, not explicitly 
written in Eq. (|2.1|) . and can be evaluated through a strict perturbative expansion (implying no 
resummations) of the full theory pressure. The second term, p^, represents the same thing for 
MQCD and can be computed through a diagrammatic expansion of the pressure of this theory, 
while the last piece, pc, corresponds to the (only) fundamentally non-perturbative contribution 
to the full theory pressure in the form of the pressure of MQCD. To order g^lng — the current 
state-of-the-art of the field — the three terms have been computed at u = in Ref. ^7] and at // 7^ 
in Ref. [q]- In Eq. I|A1|) of Appendix A we quote the result of Ref. |a| to order g^ in a simplified 
case where equal chemical potentials have been assumed for all quark flavors. 

Formally, the condition ensuring the validity of dimensional reduction is that the temperature is 
the largest fundamental energy scale in the problem, so that the above scale hierarchies are satisfied. 
In the limit of low temperatures T <C /i, the Debye mass becomes just mn = y^Nf/2gfj,/Tr which 
shows that the identification of the static degrees of freedom as the only IR sensitive ones is justified 
as long as T <C (7/U or, formally, T ~ g^n, with any x < 1. Assuming this to be the case and defining 
r = 'TTT/{g^fi), one can then extend the validity of the standard dimensional reduction results to 
the case of an arbitrary x < 1 by simply inserting this value of the temperature, 

ttT = rg^f,, (2.5) 

into Eq. (3.13) of Ref. [g] and expanding the result as a power series in g up to but not including 
order g^/j,^. The expansions of the H functions appearing in the coefficient of the 5^ term with an 
argument h/{2ttT) can be performed with the aid of the results of the Appendix D of the same 
reference. 

Here, we exhibit the result of the above procedure only for the case of x = 1/2, and as above we 
assume all quark flavors to have the same chemical potential^ /i, where we can use the simplified 



Here, one cannot simply multiply the result by -^ X] f i^i order to generalize it to an arbitrary rmmber of flavors 
with independent chemical potentials. 



result of Eq. ()A1|) . After some straightforward manipulations we obtain 
d^M'' / 4dF Sr^dp 2 f ^f lOr^ / Ydp 

^ ( 72r| In^ + 32(2Ca + STpl^V^ + 4C7ArF f 71 - 33 In 2 + 33 In - 

5767r* \ \ H 

- 153CfTf + 4r| fll - 367 + 36 In 8r - 12 In - 



q/2 2VT^{Ca+Tf)t^ 



+ 



+ ^ 3.2 

+ ^^^^(4(53C7A-6CF + 2ri.)lng + 3CF(35 + 16^^^ -161n2r 

- 227 - 20 In 2 + 487 + 104-^:) — f +1441n-^ - 1521nr - 1321n - 

C(-l) vr^ ^ 

- 4rp fs + 8 In 2 - 127 + 16 ^ , ~ , - I27 - 41nr - 121n - 

^ C(-i) /^ 

+ ^ ^ ^ - 33CAr| In^ + 4C7i7rV^ + WaTft^'^t^ + 4T|7rV^ - 54C7^r| 

12\/Tft^ V 

+ 33CyiT|,('l + 27-21n4T + 21n-) + 12r|('l + 21n2 - 21n - M 

^^^ (166CI + ^ CaTf + 72C^r^ + ^ ^^) + 0(/)|- (2.6) 

To test the validity of this expression, one may compare it with the prediction of other approaches 
entirely independent of dimensional reduction which we will do in Sec. IIVI 

B. The Preedman-McLerran result for the T = pressure 

At exactly zero temperature, important simplifications take place which resulted in the pertur- 
bative expansion of the QCD pressure being extended to order g^ in this region already remarkably 
early, in the late 1970's |9|. The most important effects stem from the fact that at T = loop inte- 
grations become much more straightforward than at finite temperature due to the Matsubara sums 
reducing to ordinary integrals. The computation of Ref. |9| was organized in much the same way 
as the one to be presented in this paper, with the exception that most of the technical details were 
much less involved. Especially with zero quark masses, the evaluation of the 2PI graphs at zero 
temperature is significantly simpler than at T 7^ 0, leaving as the real challenge the computation 
of the 'plasmon' sum, i.e. the resummation of ring diagrams. This they did in a highly imaginative 
— and to a large extent analytic — way, having to resort to numerics only in evaluating a few 
low-dimensional integrals. 

While the original calculation of Ref. 9] was performed in the momentum subtraction scheme 
and moreover using by now somewhat out-dated methods that led to largish numerical error bars, 
the result was to a large extent independently^ rederived in Ref. [^ in the ms scheme. For A^^ 



With the exception of the plasmon sum, for which similar methods as those in Ref. |9| were used. A translation 



flavors of quarks at a equal chemical potential fi, the result reads 






ATf In ^^ + - {UCa - 4:TF)ln ■ 



17 r^ 

Cf + — (568 -264 In 2) 

2 36 



- 2r^(|-141n2 + ^(ln2)2-5-^)]+0(/ln5)Y (2.7) 

where 5 ~ —0.856383209... is defined through a one-dimensional integral in Ref. 

C. Minimal HTL/HDL resummation 

At low but non-vanishing temperatures, the Freedman-McLerran result for the zero-temperature 
interaction pressure (where the Stefan-Boltzmann (SB) result has been subtracted off) to order 
g'^ needs to be corrected by T-dependent contributions. At two-loop order there are infrared safe 
terms of order g^^'^T'^ which clearly are more important than the unknown four-loop T = terms 
~ g^^i^ as long as T > g'^ii. On the other hand, when T approaches the scale g^x from the dimen- 
sional reduction region above, the leading contribution of ring diagrams Tm^ ~ g^Tfi'^ becomes 
comparable to g^fi^T^. A naive extrapolation of the Tm'^ term to temperatures parametrically 
smaller than g^ would even suggest that for T ~ g^fj, with x > 1 these become more important 
than the two-loop terms, as Tm'^ ~ g^~^^fJ. ^ g^^'^^fx ~ g^fi'^T'^. As we shall see, the T-dependent 
contributions from ring diagrams indeed become more important than the 2-loop term g'^ji^T'^ here, 
though they are not enhanced by a relative factor g~^^~^> as suggested by dimensional reduction, 
but instead only by a logarithm. 

For temperatures T < g/x, where dimensional reduction is no longer applicable, it becomes 
important to keep the nonstatic parts of the gluon self energy in the ring diagrams. At low 
momenta and frequencies of the order g^, the leading terms in the gluon self-energy are given 
by the so-called hard thermal loops (HTL) approximation with the overall m^ factor replaced by 
its zero-temperature value — a special case occasionally referred to as hard dense loops (HDL). 
In the longitudinal gluon propagator, one can observe the usual Debye screening effect at the 
frequency w <C g^, but in the transverse propagator the situation is more complicated. At strictly 
zero frequency the magnetostatic HDL propagator is massless, but for small but nonvanishing 
frequencies w <C g < nij^ its inverse has the form 

The transverse part of the propagator thus has a pole at imaginary q and \q\ = m,^^^g{uj), introducing 
a new parametrically small dynamical screening mass 



31 



/ 2 \ 1/3 

m^^^iuj) = I — - — I , a; < ruo, (2.9) 

which represents an in-medium version of Lenz's law. As soon as the temperature is small but 
nonvanishing, the ring diagrams obtain contributions involving the Bose-Einstein distribution func- 
tion which leads to sensitivity to this additional scale. In these contributions, we effectively have 



of the original result of Ref. Q to the ms scheme was previously given in the last reference of Ref. |2fil | and in 
Ref. ^. 



10 

iT^magi^ rs^ T) ~ g^'^^^'/^fi foi T ~ g^ fi and x > 1. Note that this is parametricahy smaller than 
i^D ~ gfJ-^ but always larger than the magnetic mass scale of MQCD, m„-,^^{uj = 0) = g^T ~ g'^^^/j,. 

The resummation of the nonstatic transverse gluon self-energy gives rise to terms nonanalytic in 
the temperature which to lowest order in a low-temperature expansion turn out to be of the order 
g^fjpT^lnT. This gives rise to so-called anomalous or non- Fermi-liquid behavior in the entropy 
and specific heat at low T, because instead of the usual linear behavior in T the entropy then has 
a Tin T term which is the hallmark of a breakdown of the Fermi-liquid picture (first discussed in 
the context of nonrelativistic QED by Norton, Holstein and Pincus [23]). Indeed, inspection of 
the dispersion laws of fermionic quasiparticles reveals that there is a logarithmic singularity in the 
group velocity at the therefore no longer sharply defined Fermi surface^. 

For a long time, only the multiplicative coefficient of the TlnT term in the specific heat was 
known. It was only rather recently [lO| that also the scale under the logarithm was determined 
together with the next order terms in the low-temperature (T <^ gfj) series which in addition 
involves fractional powers of T due to the cubic root in Eq. ()2.9() . For the pressure, these "anoma- 
lous" T-dependent contributions are contained in an expression, which was first derived in Ref. [la| 
and which we shall label by HDL+, 






,HDL+ _ y -^F ^Orp2 ^ / J^ „ f^ \ J.. ^2 



g 487r^ 27r3 Jq Jq 

2 „2 I TTHDL ~ 



^,, -MI!; 
2Imln ' " 



(1^-Qo 



+ Imln'^ "^0+^^^ 



+ 0{g'T^) + 0{g'i^T') + 0(<7V'T^), (2-10) 

with 5p denoting the temperature-dependent part of the interaction pressure 

(5p = Ap - Ap|t=o, 
Ap = p-psB- (2.11) 

The expression (|2.10() can be viewed as a minimal^ resummation of HDL diagrams, where the 
HDL self-energies are only kept in the infrared sensitive part of the ring diagrams involving the 
distribution function n^, while infrared safe two-loop contributions are treated in an unresummed 
form. 



1. T parametrically smaller than mo 

With 5 <C 1 and a; > 1 in T ~ g^//, the temperature is parametrically smaller than the Debye 
mass rrio ~ 5/^ and Eq. (|2.1())) contains the leading contributions to the temperature-dependent 
parts of the interaction pressure, which ignoring logarithms are of order g'^fj?T^ ~ (7^+^^|U^, while 
the higher-order terms in Eq. (|2.1()|) are at least of order g^'^'^^ ^x^. The Freedman-McLerran result 
for the T = pressure, Eq. 1)2. 7|) . is accurate to order g'^^j? and its error is of order g^ ^'^ (again 
ignoring logarithms of g). Eq. (|2.10j) thus represents the leading correction to the Freedman- 
McLerran result as long as x < 2 (i.e., T > g'^fi), whereas in quantities such as the entropy density 



A systematic calculation of the group velocity beyond the leading- log approximation has only recently been carried 

out in Ref. 0. 

As opposed to the HTL/HDL resummation considered in 1251 l26ll which aims at improving the convergence of the 

perturbative series at high temperature by retaining higher-order effects from HTL/HDL physics beyond what is 

needed from a perturbative point of view. The -I- in HBL""" and HTL"*" is meant as a reminder that the corresponding 

quantities are not expressed in terms of HTL/HDL quantities only, but combined with unresummed infrared-safe 

contributions. 
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s = dp/dT and the various specific heats, where the T = part of the pressure drops out, it is in 
fact the leading term in the interaction part for all x > 1. 

In Eq. ()2.1U|) . g appears only in the combination g'^ = g'^Tp, and it is therefore convenient to 
define a reduced temperature variable 

f = 7rr/(r/u). (2.12) 

For X > 1, the perturbative content of Eq. (|2.1()|) is that given by the low-temperature expansion 
worked out in Refs. 10|;|la|- With the above variables, this reads 

1 <5p«°^+ fV^^-i) / / 1 \ , 4 6 3\ 



Ng mi 72 



^ ^ (3) C (3) ^8/3-8(x-l)/3 , g ^ ^ (t) C(— ) -lo/3-10(x-l)/3 
3^/37r7/3 ^ 9^3^11/3 



2048 - 2567r^ - 367r* + 37r'^ _4_4(^_i) 
^ 21607r2 ^ ^ 

,14{a-l)/3^ , ^/ 2x 



In ( __^ -1 I + Invr + c 



+ 0(5^*^^-^^/^)+0(5^^'), (2-13) 

where c ~ 4.0993485 ... is given by a numerical integral defined in Ref. [la|. The latter of the 
error terms in Eq. (|2.13|) corresponds to the leading-order terms to be expected from three- and 
higher-loop contributions^ proportional to g'^fi'^T'^ which are presumably enhanced by logarithms 
of T and g. Depending on the value of x > 1, a finite number of terms in the low-T expansion 
remain more important than this (see Fig. El in Sec. IVJ) . 

When X = 1, i.e. T ~ g^, the expansion of Eq. (|2.13|) clearly breaks down (unless f <C 1) 
and the HDL-resummed expression of Eq. p.lflj) therefore needs to be evaluated numerically as in 
Ref. 16]. This expression has then the form of 5^^^ times a function of T/(gfi), and is therefore 
of the same order as the (7^ term of the T = pressure of Freedman and McLerran, to which 
it is to be added. As displayed in Ref. [la] for the case of the entropy, and as we shall see for 
the pressure in the plots of Section IIVI of the present paper, the T-dependent terms of Eq. (|2.13|1 
smoothly interpolate between a dominant g^T'^ fi^lnT behavior at low temperature and the terms 
of order g'^T'^p? , g^fi^T, and g'^fiHnT of the dimensional reduction pressure which should be the 
dominant terms at sufficiently high temperatures and which remain comparable to g'^fJ.'^ as long as 
the parametric equality T ^ gn holds. 

2. T parametrically larger than nio 

When X < 1 in T ~ g^fi, i.e. T S> ^/U, dimensional reduction provides the most accurate 
result available for the QCD pressure. Up to an error of the order of three-loop contributions 
proportional to g'^n'^T^ ~ (7^+^^/U^, one can however reproduce its prediction by extending the above 
HDL-resummed calculation to include the leading thermal corrections to the gluon self-energy. In 
practice, this means replacing the HDL approximation by the HTL one and also keeping the order 
g^T'^ terms originating from infrared-safe two-loop contributions to the pressure that were omitted 
in Eq. (|2.10 |) because they were of too high order when x > 1. This possibility was mentioned 
in Ref. [1^], but not considered further because that work concentrated on the region of T < gfi. 
For the purposes of the present paper, we however write down the straightforward extension of 



There, g^ no longer appears exclusively in combination with Tp. 
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Eq. ()2.1U|1 to the HTL approximation in the form 

J_r HDL+ ^ g^^-F 2 jn2 , g^(2C^ -Tp) ^^ 
Ng ^ 487r2'^ 288 

,2 



1 /'OO /'OO 

/ „2 _ 2 I TTHT 



2Imln 



+ n?^ 



g2 



g2 



+ 0{g''^i'T'). (2.14) 



Combining the above expression with the Freedman-McLerran result of Eq. (|2.7j) to obtain 

^pHDL+ _ ^HDL+_^^^ = Ap™^L + 5p™^^ (2.15) 

we have an expression for the interaction pressure whose error is of order (jf™™(4+2a;,6) £qj, ^t^ j, ^ g^ ^. 
This we shah compare (and thus test) in the fohowing with our new approach which resums the 
complete one-loop gluon self-energy (i.e., not only the leading HTL/HDL contribution) in ring 
diagrams. Note that the accuracy of ()2.15|) is at least of order g'^ for all parametrically small 
temperatures, excluding only the case of x = 0, where T ~ //. 

III. THE NEW APPROACH 

In this Section, we introduce our novel and strictly four-dimensional calculational scheme de- 
signed to reproduce the perturbative expansion of the QCD pressure up to and including order g^ 
at all values of ^ and T. Our guiding principle is that when faced with the necessity to sum up 
graphs with multiple self energy insertions to circumvent infrared problems, we consider the entire 
self energy and not only those parts which are identified as relevant in some effective field theory 
description, such as the Debye mass in dimensional reduction or the HTL/HDL self energy in the 
corresponding resummation schemes. Because we (at present) limit ourselves to order g^ accuracy, 
it will be sufficient to resum only one-loop self-energies in the infrared sensitive graphs, while IR 
safe diagrams will be treated perturbatively, using bare propagators. This will introduce gauge 
dependence to our results, but only at orders beyond g^ which we will explicitly discard by either 
considering values of g low enough for the higher order terms to be negligible or by performing 
numerical series expansions up to 0{g^). 

We begin our treatment with a general diagrammatic analysis where we identify the relevant 
classes of Feynman graphs that need to be considered. After that, we describe their evaluation and 
show how adding them together leads to the final result displayed in Section IIIIEI Many details of 
the calculations as well as the results of several individual pieces of the result are left to be covered 
in the Appendices. 

A. Identification of the relevant diagrams 

To determine the QCD pressure up to and including order g^ on the entire deconfined phase 
diagram of the theory, our first task is to identify all Feynman diagrams that contribute to the 
partition function at this order. These trivially include the two-^/won-irreducible (2GI) diagrams 
up to three-loop order, displayed in Fig. 1 which a straightforward power counting as well as 
the explicit calculation of Ref. [a| confirms as infrared finite for all temperatures and chemical 
potentials. 

In addition to these cases, there are, however, several other classes of IR sensitive diagrams 
that need to be resummed to infinite loop order, as a power counting reveals that the dressing 
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FIG. 1: The one-, two- and three-loop two-gluon-irreducible (2GI) graphs of QCD. The wavy hne stands 
for a gluon, the dotted Hne a ghost and the sohd hne a quark. 

of (at least some of) their gluon lines with an arbitrary number of one loop gluon polarization 
tensors does not increase their order beyond g^. These diagrams are shown in Fig. 2, where the 
first set corresponds to the well-known class of ring diagrams that leads to the known g^ and g'^ln^ 
contributions to the pressure at high T |j,y| and to the g'^lng term at T = _9]. Among others, 
this class contains the set of all three-loop two particle reducible (2PR) graphs of the theory which 
are missing from Fig. 1. 

As we shall see (in contradiction to the opposite assertion in Ref. [j|), the resummation of the 
ring diagrams is, however, not enough to obtain the entire order g^ term correctly at nonzero T. 
Although without resummation starting at orders g^, g^ and g^, respectively, the classes of Fig. 2 b- 
d, corresponding to self-energy insertions in the gluonic two-loop 2GI diagrams Id and Ig, have the 
potential to give rise to contributions of order 0{g'^T'^ ^'^) and 0{g^T'^) to the pressure. When T is 
not parametrically larger than ttt-d, it turns out that only the class b gives a non-zero contribution 
at this order, being proportional to g'^T'^m'^. When T ~ g'^n with x > 0, none of the three classes 
contributes to the pressure to order g^n'^, but in the calculation of the low-temperature entropy 
and specific heat they have to be taken into account already at order g'^fj?T. 

For any other classes of diagrams apart from those shown in Figs. 1 and 2, it is very straight- 
forward to see that the contributions will be beyond order g'^. In particular, if we were to add 
an additional gluon line with some number of self energy insertions into the graphs of Fig. 2 b-d 
(i.e. dressing the three-loop 2GI diagrams with self energies), we would notice that the two extra 
insertions of the coupling constant due to the new vertices (vertex) ensure that these graphs only 
contribute to the pressure at order g^. Similarly, one can see that the inclusion of the two-loop self 
energy into the ring diagrams only has an effect on the pressure starting at 0{g^). 

B. The 2GI diagrams 

In Feynman gauge, the sum of the 2GI diagrams in Fig. 1 at arbitrary T and /i can be directly 
extracted from from Ref. [3] with the result 



P2G1 = ^'dAT'' \^\^ + J^[l+ 30A' + 60/2^ 
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FIG. 2: Classes of IR sensitive vacuum graphs contributing to the QCD pressure at order g'^. The black dots 
represent the one-loop gluon polarization tensor given in Fig. 3a and the indices rii stand for the numbers 
of loop insertions on the respective lines. 



9(47r)2 [ "" 2 

g^ { 2?,C\ - CaTf{29 + SGO/i^ + 720 fi^) + 4r|(5 + 72p,^ + 144/2^ 



54(47r)4 1 e 

+ Ci(l821n^+247 + 272f^-90f(£|) 



_2 , ^^_4.^^ A 217 ,,C'(-1) , 72C'(-3) 



+ CaTf -16(5 + 36/i^ + 72/i*) In—— 56^:7 — f + 



47rr 5 C(-l) 5 C(-3) 

+ 24 ( 9 + 4^p^ j fi^ + 432/i^ + 144(1 + Afl'^Ml, z) + 3456H(3, z] 

+ 4T| ( (1 + 12/i2) (4(5 + 12A')ln^ + 15 + S^^^ + 36/2^ 

-2^... .,^\ n^ ^ ^35 ,,C'(-1) ,./.., -..C'(-i: 



+ 144(1 + 4^^)b^(l,z)l -9CfTpI- -16:^^ + 4(59 + 16:^^1/2^ 

+ 664/2^ + 96 (i/2(l + 4/22)H(0, z) + 2(1 + 8/2^)^^(1, z) - 12i/2H(2, z)) J I j , (3.1) 

where /2 = ^/{2ttT) and where we have renornialized the gauge coupling using the usual zero- 
temperature renormalization constant Zg. The sum, however, still contains uncanceled UV 1/e 
divergences and depends on the choice of gauge, so that it has no separate physical significance. 

C. The ring sum 

To order g^, the ring sum of Fig. 2a can be separated into three pieces pvv> Pvm and Pring 
according to Fig. 3 by decomposing the one- loop gluon polarization tensor (see Appendix [n]) into 
its vacuum (T = /i = 0) and matter parts. Note that only the matter part has to be resummed, 
as the vacuum parts contribute to order g'^ only through the two three-loop diagrams in Figs. 3b 
and cJ The evaluation of pvv and PvM is relatively straightforward, and fully analytic expressions 
for them are given in Appendix IB 11 



^ Take any graph G belonging to the ring sum and having four or more loops and at least one vacuum tensor 
insertion, and consider it in the Feynman gauge. Applying Eq. IJC4^ to it and contracting the Lorentz indices of 
the vacuum tensor with one of its neighboring gluon propagators, we see that G is proportional to g^ times a similar 
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FIG. 3: a) The one- loop gluoii polarization tensor n^y(P) divided into its vacuum (T = /i = 0) and matter 

(vacuum-subtracted) parts. 

b) The IR-safe Vac- Vac diagram contributing to the pressure at 0{g'^). 

d) The IR-safe Vac- Mat diagram contributing to the pressure at O(g^). 

d) The remaining 'matter' ring sum. 

To evaluate the remaining matter ring sum p^i^^ we define the standard longitudinal and trans- 
verse parts of the vacuum-subtracted polarization tensor at d = 4 — 2e by 



u,{p)6^'' = ^(ng[i(p)-ngg(p) 



p. 



n^{p)5' 



ab 



d 



KiP) - K(P) ivac -^ «iP) - KoiP) 



p2 

p2 



vac / I ) 



(3.2) 
(3.3) 



where we have used the fact ^363 that the one-loop gluon polarization tensor is transverse with 
respect to the four- momentum P in the Feynman gauge. In terms of IIt and IIl, the sum of the 
ring diagrams is then readily performed with the result 



Pr: 



■H 



In 



l+U^iP)/P^ -U^iP)/P' 



+ (d-2)(ln[l+nT(P)/p2] - U^{P)/P^' 

which is now explicitly IR safe. 

As the functions nL(P) and Ilrr{P) behave at large P^ hke (see Sec. B.1.2 of Ref. |l7| ) 



(3.4) 



n 



L/T 



{P) ~ > -2{l + e)CA9'i^+0{l/P') ^ n„v + 0(l/P'), 



(3.5) 



it is, however, immediately obvious that the sum-integral of Eq. (|3.4j) is still logarithmically di- 
vergent in the ultraviolet at T 7^ 0. To regulate the divergence, we add and subtract a term of 
the form {\ + d — 2)(nuv)^/(2(-P'^ + m?)'^) from the integrand, with m being an arbitrary mass 
parameter shielding it from IR divergences. By further adding and subtracting the corresponding 
massless term from the counterterm, we obtain three separate contributions to Pring^ an UV and 
IR finite (at least to order g'^ — see below), ?TT--dependent ring sum p^^^l", an UV finite, but IR 
divergent and m-dependent p™g and an UV and IR divergent and massless p^^^ 

= „fl"it<= + „™ +„uv 



graph with the vacuum insertion removed. But this graph is nothing but one of those diagrams that appeared in 
the original sum which implies that G has to be proportional to at least the fifth power of the coupling. 
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^.^^^^¥^i\^.A^-^\ (3.7) 



Hnito ^ (^A 

2 

+ 2 fin 



48 Tp\ (P2 + m2)2 pi 



■faring 4 



i(d-l)(iA(nuv)'^^. (3.8) 



The two first terms can be evaluated numerically at e = while the third one needs to be regulated 
with finite e. It is noteworthy that one can set e = even in the formally divergent p™g due to the 
fact that its IR divergence originates solely from the zeroth Matsubara mode of its second term 
which vanishes identically in dimensional regularization. The explicit values of plf^^ and p^^^ are 
given in Appendix IB 21 while the numerical evaluation of p^^^^" is the subject of Appendix IdI 

D. The double and triple sums 

If the sums in Figs. 2b-d were to start from n = 0, these multiple resummations would clearly 
correspond to the dressing of the propagators in three two-loop diagrams with the one-loop gluon 
polarization tensor. In the present case, we instead define a four-dimensionally transverse^ tensor 
/S.G^u{P) by the equations 

^^"■^^^ " p2 + nL(P) " p2 = -p2(p2 + nL(P)) ' ^^'^^ 

^^"^^^ = p2 + n,(p) " p^ = -p2(p2 + n,(p)) ' ^^-^"^ 

corresponding to the difference of a dressed (with the vacuum-subtracted self energy) and a bare 
gluon propagator in the Feynman gauge. It is a straightforward exercise in combinatorics to show 
that the symmetry factors of all graphs in Figs. 2b-d equal 1/4 independently of n — a result 
particularly obvious in 2PI formalism. To order g^, these three classes of diagrams, denoted here 
by Pbi Pc and p^, can then be written in the forms 

dACA 24 ^G^^,{P)^Gpp>{Q) 



^' 4 ^ Ji,Q {P + Qf 

X {g^'^ilP + Q)P - g''P{2Q + Py + g^^iQ - P)' 

X (g^^'-{2P + Q)p' -g^p'{2Q + PY' +gp'^''{Q-PY]+0{g^), (3.11) 



^^ TpQ 

X {g^"'{2P + Q)P - g''P{2Q + Py + g^^iQ - PY) 

X f5MV'(2P + Q)p'_/p'(2Q + P)M'+/M'(g_P)-')+0(/), (3.12) 



P, = -^^g^i (AG^^(P)AG,,(Q)-AG^,(P)AG^,(g)) + 0(/). (3.13) 

^ Tpq 



Thus decomposable into three-dimensionally transverse and longitudinal parts. 
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All contributions involving the vacuum piece of the polarization tensor have been discarded as 
being of order g^ , following a reasoning similar to that in Footnote [3 

It is worthwhile to first perform a power counting analysis to determine at which order the 
above sum-integrals start to contribute to the pressure. In the regime of dimensional reduction, 
where T ~ g^ ^ with x < 1, one merely needs to consider the contributions of the zeroth Matsubara 
modes, as for the others the temperature acts as an infrared cutoff, leading to their values being 
proportional to at least g^Tfi^ ~ g^'^^n'^. In the region T ~ gn, the Debye mass is, however, 
of the same order as the temperature, implying that all Matsubara modes give contributions to 
the pressure parametrically similar in magnitude. Scaling the three momenta in the integrals of 
Eqs. H3.11() - ()3.13() by gfj,, one quickly sees that the results for the sum-integrals in this regime can 
up to 0(5^) be written in the form g'^T'^^^ f{T/{g^)), where the contributions of the non-static 
modes to the function / vanish as the parameter T /{g^) approaches infinity, while in the opposite 
limit T/(gfi) -^ the function approaches a constant. As long as we are interested in the value of 
the pressure only to order g^fi'^, these graphs can clearly be altogether ignored. They will become 
relevant in the determination of the 0{g^ii^T) contributions to the specific heats, but this is outside 
the scope of the present work. 

For now, we can concentrate our attention to the regime of dimensional reduction and therefore 
to the zero Matsubara mode parts of the above sum-integrals. Here, we encounter an important 
simplification which results from the fact that only the longitudinal part of the static gluon po- 
larization tensor has a non-zero zero momentum limit at one-loop order. As the finite momentum 
corrections to the functions nL(P) and nx(-P) clearly correspond to higher perturbative orders, we 
can simply replace 

2 

Tfl 

^GAP) - - p2(p2;^2) '^MO<^.o (3.14) 

in the integrals, leading to a dramatic reduction: both p^ and p^ then vanish identically. This 
can, however, be easily understood from the point of view of the three-dimensional effective theory 
EQCD as a demonstration of the fact that the ^q and Aq operators in its Lagrangian are not 
accompanied by couplings of order g and g"^, respectively, but only g^ (at nonzero /i) and 5^. 

In contrast to the above, for p^ one does obtain a non-zero value which has a direct parallel in 
EQCD in the form of an 0{g) coupling between one massless Ai and two massive Aq fields and a 
corresponding two-loop diagram with one Ai and two ^0 lines. Applying the limit of Eq. (|3.14|) to 
the sum-integral of Eq. (|3.11|) , it is easy to see that we can reduce the expression of p^ (to order 
g^) to the simple form 

^ dACA 2^, 2 fj^ [^ (p-q)' .3l5^ 

^' 4 °^ y (27r)3y (27r)3p2(p2 + m2)q2(q2 + m2)(p + q)2 ^"^'"^^ 

which can be solved straightforwardly by introducing three Feynman parameters and using stan- 
dard formulae for one-loop integrals in three dimensions. After some work, we get 

dACA^2 4 2/ f d?p 1 f' f d'q 1 



(27r)3 p2(p2 -\- m^) Jq J (27r)3 q2 -|- 2xq • p -|- xp2 -|- (1 — x)m'^ 



\ f <fp 1 

ax I 







(27r)3 (p2 + xm2 )2 



d.c. rV„.^ fi t\. 1 f-. ' , ' /'■ > 



4 (47rmD)2 [vr Vo yJx{l-x)J\i y/V 'i- - V + v/x Jq ^ 
dACA ^2_2 g^ 

(4^) 



T^m^^-^ (1-4 In 2) (3.16) 
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which we identify as the entire contribution of the classes b-d of Fig. 2 to the QCD pressure up to 
order g^. 

E. The result 

We are now ready to write down our final result for the pressure, valid on the entire deconfined 
phase of QCD and accurate up to and including order g'^. Assembling all the various pieces, this 
function reads 

P = ip2Gi + Pvv + PvM + P^, + P.) + (pZ, + Pfinr) + <^i9'Tf^') + 0(gV4) (3.17) 

= Panl+P:rf, + 0(/V)+0(<7V'), (3.18) 

where p^^i stands for the sum of the first five terms in Eq. 1)3. 17|) and 

p:-:; = p^::t+p1, (3.19) 

is to be evaluated numerically. One should note that in this notation all 771-dependence in contained 
in the two pieces of p^f^g; naturally canceling in their sum. In addition, it is worthwhile to point 
out that the inclusion of the term ph in Eq. p.l7|l is inconsistent in the region of T ~ g^n, x > 1 
where we have neglected several contributions of the same magnitude. As this term, however, is of 
order g^^'^'^n^, i.e. at least of order g^fj.^ in the region in question, the inconsistency is in any case 
beyond the order to which our result is indicated to be valid and can therefore be ignored. 

Collecting the expressions for all of its parts from above and from Appendix ^ the function 
p,„i reads 



n4/ 1 f-, , dp M 






')] 
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27(47r)4 [ ^V 47rT ' ({-1) C(-3) 



-2 , -,.o.-4M A , 2391 , r'(-l) , 116C'(-3) 



C.T.^(47 + 792,^ + 1584,^)ln^+^+4|^+ 5 C(-3) 

,C'(-1)A.2 , ooon-4 , .no.-, , . -2^ 



+ 6 ( 257 + 88 ^' ' J /2^ + 2220/i^ + 792(1 + 4/Z^)H(l, z) + 3168H(3, z) 
+ T^( (1 + 12P^) (4(5 + 12,= ).nA + W^^) + ^ + '' ^''"'' 



attt c(-i)y 55 c(-3) 

+ 312/2^ + 624/2^ + 288(1 + 4/2^)^(1, z) + 1152H(3, z) 



+ - CpTphb - 32(1 - 4/2^)^7 — r + 472/2^ + 1328/2^ 

+ 192 (i/2(l + 4/2^)N(0, z) + 2(1 + 8/2^)^^(1, z) - 12z/2b^(2, -z)) ) I ) • (3-20) 

Not only have all the UV divergences canceled between the different parts of this result, once 
the renormalization of the gauge coupling g has been taken care of, but this expression actually 
contains all the (explicit) renormalization scale dependence of the pressure up to the present order 
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in perfect agreement with Ref. |a|, leaving pl^^"^ entirely independent of the parameter A. Eq. ()3.2U|1 
is also valid for all values of T and fi; the limit for // — > is given in Eq. (|B7|) and the limit 
T — > in Eq. HB8|) . All terms non-analytic in g'^ are contained in the piece p^f^" awaiting numerical 
evaluation. 

In the following we shall denote our final result for the pressure — which is accurate to order 
g^ for all values of T and /i (while also containing some incomplete contributions of higher order, 
to be discarded later) — by 

Piv = p.„i+p;r.:;- (3.21) 

F. Numerical infrared issues 

Before moving on to examining our result by numerically evaluating the function jp!^"'^" in 
Eq. 1)3. 6|) . there is one more practical issue related to the magnetic mass problem |28l . |29| that 
needs to be dealt with. To wit, in the limit P ^ 0, the argument of ln(l + IIt^[P) / P'^) becomes 
negative, resulting in an unwanted imaginary contribution to the integral which actually renders 
Pring° infrared singular beyond order g'^. This problem depends on the choice of gauge, but is 
present in all covariant gauges (as well as the Coulomb gauge). 

The origin of the problem can be traced back to the fact that when dressed with the full one- 
loop self-energy, the transverse part of t he g luon propagator develops a space-like pole. For po = 
this pole is determined by the equation [22| 

p2 + n,(po = 0,p) = p^- g^N,T^±^ip^p (3.22) 

where £, is the gauge parameter of covariant gauges.^ It is evidently unphysical and appears only 
at the non-perturbative magnetic mass scale g^T, which contributes to the pressure starting at 
order g^T^. This suggests that we can in fact eliminate the entire problem by adding by hand a 
magnetic mass term to the transverse self-energy in Eq. (|3.6|) 

n,(P) ^ n,(P)+mL, (3.23) 

with (for ^ = 1) 

m^^, = Cf^g^CAT (3.24) 

and Cf > 1, which only has an effect on the pressure beyond 0{g^). Indeed, comparing with the 
effective magnetic mass for nonzero frequencies, Eq. ()2.9() . we find that the magnetic screening 
behaviour is modified only for frequencies po ^ 9^T when /i ~ T and even pq < g^T{T'^ / fj?) when 
T -^ II. Note, however, that the introduction of this magnetic mass for the transverse self-energy 
alters the UV behavior of pj^'^g", implying that both p|^°g° and p^^^ have to be modified to account 
for this reorganization. In pf^g", this change is crucial because it renders the result finite, but for 
the already finite p|;f^g the effects are beyond the order of interest (see Add. IdI) . 

The numerical evaluation of p^^^l" is performed along the lines of Refs. |15l . l32| . with the sum 
over Matsubara frequencies being converted to an integration in the usual way (see e.g. Ref. 33]). 
Contributions containing the bosonic distribution function ref, are best evaluated in Minkowski 



^ Replacing the ordinary one-loop gluon self-energy by one that includes resummation of the Debye mass does not 
cure the problem, but only produces a different gauge-dependent spacelike pole |3(ll3ll|. 
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space, as UV problems are cut off by rij,, while the other contributions are evaluated in Euclidean 
space in order to numerically exploit the Euclidean invariance of UV contributions. By varying 
the parameter Cf in Eq. ()3.24|) , we can verify that the effects of this infrared regulator are indeed 
beyond the order g we are aiming at. The remaining part, however, gives rise to yet another type of 
unphysical pole, which (at least in the long-wavelength limit) has been well-known since the earliest 
perturbative calculations in finite-temperature QCD |2J]: in covariant gauges, the one-loop gluon 
self-energy, evaluated at the location of the poles corresponding to time-like propagating plasmon 
modes, gives rise to a (gauge-dependent) damping constant oc g^T with negative sign (for all gauge 
parameters ^, though not in Coulomb or axial gauges J35l. l3fi| ) . A consistent systematic calculation 
of the plasmon damping constant to order g^T requires the use of a HTL-resummed gluon self- 
energy which finally leads to a positive and gauge-independent result jsil, |22|- The corresponding 
pole is then on the unphysical sheet where it would cause no problem for the evaluation of p^^^"- 
With the bare one-loop gluon self-energy appearing in our integrand we, however, have poles on 
the physical sheet, connected to the light-cone by a branch cut, and we need to avoid them by 
deforming the contour of the numerical integration in Minkowski space as sketched in Fig. EJ The 
details of this procedure and the entire numerical calculation are described further in AppendixlDl 

IV. NUMERICAL RESULTS 

Having the result of Eq. (|3.17|) for the QCD pressure now finally at hand, we move on to examine 
it numerically by evaluating the function p^f^J" using methods reviewed in Appendix ^ and adding 
to it the analytic part of Eq. (|3.2flj) . The sum total we call piy as a reminder that its accuracy is 
of order g^ for all T and /u, while it also includes incomplete and gauge dependent higher-order 
contributions. For the most part of the following analysis, we shall explicitly eliminate the latter 
effects by either considering sufficiently small values of g or performing numerical expansions of 
our results in powers of g. 

We begin by inspecting the region where the temperature is parametrically larger than the Debye 
scale and the results of dimensional reduction should be applicable, then continue towards making 
contact with HDL results on non- Fermi liquid behavior at T < ttt-d, and finally the Freedman- 
McLerran result for T — > 0. In all plots of the present section we use the values N^ = 3, N{ = 2. 
Because of the latter, we conveniently have Tp = 1 and therefore t = f for the reduced temperature 
variables introduced in Eqs. (|2.5() and ()2.12() . respectively. 

A. T parametrically larger than ttiq 

The first non-trivial check on our result — and that of dimensional reduction — is to verify that 
their predictions for the pressure agree to order g'^ for all temperatures and chemical potentials 
that are of equal parametric order in g. This is particularly important in order to clarify that the 
(entirely correct) statement in the literature about dimensional reduction being valid as long as 
ttT is the largest dynamical energy scale does not imply a condition ttT > fi, but rather ttT » rrio 
(or even ttT > ttt-d, as we shall find to be sufficient below). To this end, we start from the most 
widely studied region of ^ = by comparing our numerical result to that of the analytic one of 
dimensional reduction, and then increase the 0{g^) value of fi/T up to /^ S> vrT while still having 
ttT » rriD ~ gfJ.. 

The results of this comparison are shown in Fig. 01 where we plot the order gHng and g^ 
contributions of the ring sum of Eq. 1)3. 6() to the pressure together with the same quantity extracted 
from the result of dimensional reduction (obtained by subtracting the analytic part of our result 
from the DR one). The agreement is perfect up to the numerical accuracy of our result, and only 
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FIG. 4: Comparison of the g^ log g and 5"* terms of the numerical computation and the analytic DR result, 
for various values of /i/T. The perturbative terms are subtracted up to order g'^ . 
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FIG. 5: Comparison of the HTL+ pressure and our numerical result piv in the region of T = r^fi, r = 0.2, 
with the known perturbative terms from dimensional reduction subtracted and the entire quantities divided 
by g^'^. This plot shows that both the HTL-t- result and our numerical one are accurate at least up to order 
5^'^. The renormalization scale has been varied between /i and 4/i. While piv — Pdr is scale independent, 
Phtl+ — Pdr has a scale dependence at order g^fi^T^ ~ .9^/^^- 



at larger values of g can one see that the agreement is getting slightly worse with increasing /x/T. 
This was, however, to be expected, since there n/T is coming closer to the value g~^ , making 
rriu/T of order one which is parametrically the limit of applicability of dimensional reduction. Our 
conclusion is that the result of dimensional reduction is valid at in principle arbitrarily large 0{g^) 
values of ///T, though the expansion in g only makes sense at smaller and smaller values of g as 
this parameter is increased. This statement will be made more concrete in the following sections. 
The logical next step is to test the validity of dimensional reduction at temperatures larger than 
but now parametrically closer to the Debye scale. For concreteness, we specialize to the case of 
T ~ y^/U, for which the prediction of dimensional reduction is given in Eq. (|2.6() . In this region. 
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FIG. 6: Same as Fig.jSJ but normalized to g^ . While the HTL+ result is no longer accurate to this order 
and diverges logarithmically, our numerical result still correctly reproduces the dimensional reduction result 
for the pressure at order 5^^'*. 



the error in our result is of order g^^'"^ fi^ and that of the minimal HTL resummation (7^/i^, so that 
the first one should be able to reproduce the first seven and the latter the first six terms of the 
series (|2.6() . And indeed, a numerical evaluation of both Eqs. H3.17() and H2.14I) and the subtraction 
of the first terms of Eq. (|2.6j) shows the expected results: as displayed in Fig. [3 we find perfect 
agreement in comparing the dimensional reduction result with the HTL one H2.15() and with that 
of our new approach up to order g^'"^ . In Fig. |HJ we see that our numerical evaluation of piv is 
accurate enough to even verify the g^^i^ term in the dimensional reduction result, while the HTL 
result starts deviating from the DR one at this order. 



B. T comparable to ttid 



In Figs. [THHl we plot the the temperature-dependent contributions to the interaction pressure 
6p (see Eq. (|2.11j) ) for T ~ ttz^ as extracted from our numerical calculation of piv but with no 
expansions in powers of g. We compare this with Phtl+ as well as with the dimensional reduction 
result expanded to orders 5^, g^ , g^ and g^ which refer to the counting in powers of g when 
T ~ /x. For T ~ gii, however, the terms g'^fi'^T'^, g'^fi^T, and g'^nHnT all become of the same 
order of magnitude and together constitute the leading temperature-dependent contribution to 
the interaction pressure p — psB which is contained in the result marked by the dashed line "g'^" . 
For completeness, we also include the complete dimensional reduction result to (explicit) order 
g^, but it should be remembered that the term g^Tfj.^ is already of the same magnitude as the 
unknown g^fi'^ piece when T ~ gfi, and is therefore both incomplete and beyond our scope which 
also explains why the g'^ curve seems to produce better agreement with our results than the g^ 
one. 

The different results are normalized to the leading term of the T-dependent part of the inter- 
action pressure in the dimensional reduction result HA1|) . 
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To understand the structure of these figures, note that the g^ curve goes like — 1 + (4/37r)?TT-D/T 
for small T and like —1 -|- 1.075 foi' large T. At T <C rriu it, of course, deviates from the exact 
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FIG. 7: Thermal contribution to the interaction pressure 5p as a function of T/m^^^ for fixed chemical 
potential /^ and coupling g = 0.1. For this value of the coupling, the results of the numerical evaluation of 
Pani + Prfng ^-^id HTL^ coluclde within plot resolution. The result is compared to the dimensional reduction 
pressure at orders g^, g^ , g'^, and g^ (where the latter is included only for completeness, as neither piv nor 



Pjjrp]^+ contain contributions of order g^) 
not visible for this value of the coupling. 



The effect of varying the renormalization scale Aj^ = ^ ... 4/i is 
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FIG. 8: Same as Fig. |3 but for g — 0.5. The results of the numerical evaluation of Pani + Pltns ^^^ htl+ 
can now be distinguished due to their different content of higher-order terms. When two lines of the same 
type run close to each other, they differ by changing the renormalization scale A^^ — fj, ... 4/i. 



result which is instead dominated by the leading g lnT~^ behavior of the low-temperature series 
of Eq. (|2.13|) when normalized by the absolute value of Eq. H4.1() . 

For small values of g ~ 0.1, Fig. shows that the numerical evaluation of Eq. (|3.17j) perfectly 
agrees with the result of the HTL resummation (the two curves lie virtually on top of each other) . 
At this value of g, also the complete dimensional reduction result to (explicit) order g^ is virtually 
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FIG. 9: Same as Fig. |S1 but with piy separated into Pani and p^fn^. As the g'^ contribution in (5pani only 
amounts to a small correction (of effective order g^), the shape of the full pressure curve as a function 
of T (beyond the rather trivial g^ contribution) is mainly determined by p^fig- The renormalization scale 
dependence Aj^ = /i ... 4/i is entirely due to Pani- 



indistinguishable from the order g^ result. The dimensional reduction result reproduces the numer- 
ical results remarkably well down to temperatures of about 0.2 rnf^^^, but at even lower T severely 
overestimates the logarithmic growth of dp/T"^ as T — > 0. This is to be expected, since, in the limit 
T — > 0, the plasmon term of order g'^n^T in the pressure is clearly unphysical, as it would lead 
to a nonvanishing entropy at T = 0; the g'^fj.HnT term of the dimensional reduction result, while 
evidently crucial for good agreement down to T ~ 0.2mD, would even lead to a diverging entropy 
as T — > 0. The point at which the dimensional reduction result ceases to be a good approximation 
for both Phtl+ and piy seems to agree rather well with the value of T/niYy where 6p switches sign. 

In Fig. ISlwe consider a larger coupling g = 0.5, for which we begin to see effects from varying 
the renormalization scale A^^ in our result by a factor of 2 around the central value 2^, except 
in the HTL-|- result, where Ky^ appears only in the T = (Freedman-McLerran) part of the 
result.^'' For small T/rn^^^, we find good agreement between the HTL+ result and piy, with the 
dimensional reduction result to order g^ lying in between the two in the range T/rn^^^ ~ 0.1 . . . 10, 
but deviating again abruptly for T/m^^^ < 0.2, which is where 6p changes sign. At this value 
of the coupling, the complete order g^ result of dimensional reduction is still reasonably close to 
the order g^ result. While it is certainly unreliable when 6p > 0, the order g^ result suggests that 
taking into account the next higher orders in g may move the onset of non- Fermi-liquid behavior 
to slightly larger T /m-Q. 

Fig. 121 shows how the final result Spis, is composed of the infrared-safe piece Pani and the ring 
sum pjifng- At parametrically small T ~ g^, the T-dependent terms in the interaction part of Pani 
which are of effective order g'^^'^ come just from the terms g'^T'^iJ?, so that the shape of 6p in the 



^^ In Fig. |H] the value g — 0.5 is kept fixed for all A-^jg, which means that the a;-axis does not correspond to a 
renormalization-group invariant variable. The (explicit) dependence of the results on Aj^ is here shown only to 
assess the theoretical error in the numerical comparison between the different approaches. Taking into account 
the implicit Aj^ dependence of g, the scale dependence of all the results we are comparing is of the order of their 
error, which is 0{g^fj^) at T ~ g/i. 
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FIG. 10: Same as Fig. [7| but for 5 = 1. At this value of the couphng, the numerical result for piv begins 
to be visibly affected by the choice of the magnetic mass (|3.24() which here is taken with Cf — I. 



(2) 

above figures is mainly determined by 5p^f,J° which is seen to coincide with (5phtl+ — <^Pdr up to 
terms beyond g^ accuracy. 

Finally, in Fig. ^J we consider g = 1 which is roughly the value of the QCD coupling at 100 GeV. 
Here the result for pj^ still follows Phtl+ for T < ruj-, and the g'^ result of dimensional reduction 
for T ^ ?TT,D, but there is an overall shift due to the (incomplete and gauge dependent) order g^fi'^ 
terms in p^y. At this value of g one also notices that the dimensional reduction result to order g'^ 
deviates rather strongly from the result to order g^, thus showing poor apparent convergence, in 
particular at high r. 

For the remainder of the discussion of our numerical results at T ~ mo , we concentrate on the 
contributions of order g^fi'^ to the pressure and explicitly discard all terms beyond this accuracy, 
as this helps us to better analyze the breakdown of dimensional reduction observed in the three 
previous plots. To this end, we consider the difference of the g'^p-'^ term in the DR result at 
T ^ gfi and the corresponding piece in the Freedman-McLerran zero-temperature expression. A 
straightforward evaluation of this quantity gives 
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7687r6 



[33 -3d- 127 - 27r2 + 2 In 2(8 In 2 + 7) + 12 ln(7rf )] + 0{g^ 



which, like the HDL+ result, only depends on g and T to order g'^ through the combinations 
g = Tp g and r = TTT/{gp) = T/m-o and where 5 is the numerical constant appearing in Eq. ()2.7() . 
As we have seen, in this region the g'^ content of our pi^ agrees perfectly with that of Phdl+ defined 
by Eqs. ()2.10|) and (|2.11|) . Therefore, to simplify our numerical efforts, we compare in Fig. ^J 
the above function with the HDL"*" result for the same quantity 5p^ after by hand subtracting the 
terms proportional to r^ and r of Eq. (|4.2() from both results (which correspond to the terms of 
order g'^p^T'^ and g^iJ?T in the expansion of the pressure). ^^ This helps us to expose the term of 



Note that these terms are not present in the HBL"*" result at small values of t (but that this does not matter, as 
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FIG. 11: Plot of the last term in Eq. (|4.2|) (dashed line) in comparison with (5phdl+ with the first two terms 



of Eq. (|4.2() subtracted, in units of ; 



g /i /tt and as a function of log]^Q(r). 



order g^fiHn T, whose divergence at T — > signals the failure of dimensional reduction to correctly 
describe the zero-temperature limit of the pressure. 

From Fig. llll we observe that for r > 0.2mD there is a In r term in the HDL"*" result which agrees 
perfectly with that of Eq. ()4.2|) . This fact turns out to have a natural explanation which gives us 
important insight into the breakdown of dimensional reduction. The key observation is that the 
ultimate reason for this breakdown lies in the incorrect treatment of low-temperature IR divergences 
in DR: in deriving its prediction for the QCD pressure, one assumes that the temperature acts as the 
sole IR cut-off for all non-zero bosonic Matsubara frequencies in the logarithmically IR divergent 
three- loop ring diagrams, so that for them no resummations are necessary. ^^ While this indeed 
is justified for T ^ rn-n, there is an obvious problem in the region T ~ rriu, ultimately leading 
to the diverging of the DR result in the zero-temperature limit where discrete Matsubara modes 
no longer exist. In a physically consistent calculation, where one performs a resummation also 
for the non-static modes, the logarithm of temperature in the dimensional reduction result above 
gets replaced by one of a true IR regulator T x R(gfj,/T), where the function R is linear at large 
values of its argument and approaches a non-zero constant as gfJ,/T -^ 0. Unlike T, this regulator 
therefore does not vanish at T = which explains the finite r — > limit of the HDL"*" curve in 
Fig. 1111 However, at large values of f the logarithm of TR gives rise to the Inf behavior with the 
same coefficient as in the dimensional reduction result which is what we observe in Fig. IllL 

C. T smaller than mu 

In the limit f = T/m-a -^ 0, our evaluation of piv approaches the Freedman-McLerran result 
1)2. 7|) as shown in Figs. ^1 and EJ where the difference of Piy{T = 0) and Pfmcl is plotted as 
a function of g. The agreement to order g'^ and the absence of a g^ term is shown in Fig. I12| 
where we normalize the result by g^fj,^. Fig. ^Jon the other hand shows that piy differs from the 



their effect in Fig. 11 in any case vanishes as r ^ 0). 
^^ Which amounts to expanding the already resummed propagators for the non-static gluon modes in these diagrams 
in powers of gji/T. 
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FIG. 12: Comparison of our numerical result pw and the Freedman-McLerran result. The plot shows 
agreement to order g^ , i.e., agreement in the nonvanishing coefficients up to and including order g*, and the 
absence of order g^ contributions. 
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FIG. 13: Same as Fig. ^1 but divided by g^ instead of g^ , revealing that piv contains a term of order 
g^lny/i"*, which is however incomplete as it is beyond the accuracy of our setup. 



Freedman-McLerran result at order ^r^, where it contains a term of order g^\n.g. This, however, is 
incomplete, as the true (unknown) g^\ng term in the T = pressure gets contributions also from 
the two-loop gluon self energy. 

In Fig. E] we plot the coefficients of the low-temperature pressure in the perturbative expansion 
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as a function of r = TrT/{gfi) for Nc = 3 and Nf = 2. Here, we have again numerically evaluated 
the HDL+ result which as we have shown agrees with our new approach to order g'^ for T ~ gfi. For 
r = the coefficients reproduce the Freedman-McLerran result, whereas for larger r the coefficients 
are dominated by Stefan-Boltzmann contributions T^/i^ ~ 5^^"^ and T^ ~ g^fi'^. 
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FIG. 14: Dimcnsionless coefficients of tlie expansion of the pressure in powers of g at ttT = rgji as a 
function of r for Nc = 3, Nf = 2. The coefficients po and p'^ are constants, p2 is quadratic in r, while p^ 
shows a T^ logr behavior for smaU t. The renormahzation scale is varied through Ajjjg — fi . . . 4/j,. 



In Fig. 1151 we compare the T-dependent contribution to the interaction pressure, 5p, with the 
extrapolation of the dimensional reduction result (where Ppmcl has been subtracted since pun does 
not exist at T = 0) to T ^ tti-d in a linear plot. Here, 5p is normalized by m^ ~ fi'^A*^ rather 
than divided by a term quadratic in temperature to better show its absolute magnitude^'^. As 
could already be seen in the previous plots, the dimensional reduction result works surprisingly 
well down to temperatures where 6p changes sign which happens at T ~ 0.2mD. At even lower 
temperatures, the dimensional reduction result overestimates the leading positive contribution to 
Sp which in reality goes to zero like (or grows ~ logT with the normalization of Figs.EHSl) 



6p{fi, T <C gfi) 



^,VT^ln^ + 0((,^ 



(4.4) 



with the constant c ~ 0.284794. This expression determines the total pressure of the normal QCD 
phase as 



Pnormal(^, T < fir/i) = Pnormal(/X, 0) + 5p{fl, T < g/i) , 



(4.5) 



where Pnormai(Ai,0) is given in Eq. (|T7|) . 

For exponentially small T with ln{gfi/T) > g~^ the contribution of Eq. 1)4.4(1 can even be larger 
than the leading T-dependent^^ term oc fi^T^ in the Stefan-Boltzmann pressure |38l |. However, be- 
fore one reaches temperatures so small that ln{gn/T) > g~'^ in a nonabelian plasma, one encoun- 
ters the nonperturbative pair instability of color superconductivity and the formation of a gap [lj| 
~ Tc ~ ;ug~^ exp[— 37r^/((7-v/2)] at the parametrically larger temperatures with ln{g^/Tc) > g~^. 
In our (resummed) perturbative approach, we will not directly encounter this nonperturbative 
instability, but since at the superconducting transition temperature the pressures of the normal 
and superconducting phases in any case have to be equal, it is of some interest to evaluate the 



^^ The normalization to Spj^^^ is however more appropriate if one is interested in the magnitude of the effects to 

entropy and specific heat. 
^^ Note that in the entropy or in the specific heat these terms constitute the leading ones. 



29 



5p/{dA{m^ 



T=0\4 



0.002 



0.001- \ 



-0.001 



-0.002 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 

r=o 



T/ 



m 



FIG. 15: Plot of the T-dependent part of the interaction pressure 5p (see Eq. (|2.1Hl 'l to order g'^ in the 
regime ttT < TOd in units of (tti^^'^)'* = g^ii'^/n^ and as a function of f = T/mu- The dashed hne denotes 
the dimensional reduction result to order g*, and the full line the HDL+ result which in this regime coincides 
with the order-g^ content of piv- 



correction of Eq. H4.4|) at the transition temperature T^*^. Inserting ln{T/fi) ~ — 37r^/((7v2) into 
Eq. (|0|) . one finds 



PnormaU/^) -ic j " Pnormal (,/^) " j + A^ ^ ^ ^ C J ' 



(4.6) 



Comparing this to the T = pressure of, e.g., the 2SC phase from [39| and expressing it in terms 
of Tc by replacing (p — > TTe~'^T^'^, we have 



Psc(/", 0) ss p^^,^^i{n, 0) + e "^^TplJ^iT, 



SC^2 



(4.7) 



Thus we observe that at the temperature where color superconductivity sets in, the contribution 
of resummed perturbation theory is gdA^^'^ /{'^^V^) ^ 0.75(7 times the T = contribution of the 
gap. 

V. SUMMARY OF PERTURBATION THEORY ON THE /x-T PLANE 



As we have seen in the previous Section, the weak couphng expansion of the QCD pressure goes 
through changes in its form when T/ ^ becomes comparable to some positive power of the coupling 
constant g and this power is then increased. In particular, at T < g^ dimensional reduction 
ceases to be applicable and a resummation of the nonstatic parts of the gluon self energy becomes 
necessary, although numerically the dimensional reduction result works surprisingly well down to 
the rather small value T/m-Q ~ 0.2. At this temperature, which to leading order in the coupling 
reads 



Tnfl 



0.064^^5, 



(5.1) 



the T-dependent contributions to the interaction pressure change sign, marking the onset of so- 
called non-Fermi-liquid behavior. At parametrically even smaller temperatures, one eventually 
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FIG. 16: The dividing line between the regime of dimensional reduction (DR) and that of non-Fermi-liquid 
behavior (NFL) as given by Eq. (|5.1|) for iV/ — 3 (full lines) and Nf = 2 (dashed lines), in comparison 
with the weak-coupling result (|5.2|) for the critical temperature of color superconductivity (CSC) when 
extrapolated to large coupling. 



encounters the critical temperature of color superconductivity, which has been calculated to leading 
order in a weak coupling analysis as [40| 
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for a spin-zero condensate (which gives the largest value of Tc). 

In Fig. El we compare Tnfl and T^''^ when extrapolated to large coupling (where of course 
correspondingly large modifications can be expected). At least for smaller coupling, i.e., sufficiently 
high densities, there is a clear separation of regimes in the ^-T plane with qualitatively different 
weak coupling descriptions. 

In Fig. El we display the structure of the weak coupling expansion by showing how the mag- 
nitudes of its terms depend on the power x in the order-of-magnitude equality T ~ g^fi (which 
of course removes the superconducting phase from the picture). The x-axis in this figure corre- 
sponds to the ratio ln(T / fi) /ln{g) in the limit g ^ and T/fi -^ 0. For each value of x we give 
the orders of the first several terms in the expansion of the pressure in powers of g, not counting 
separately terms with an extra factor of lm.{g). Full lines denote known contributions, while dashed 
and dash-dotted lines correspond to the as-yet unknown ones. 

At X = 0, i.e. T ~ /i the situation is still the same as with /i = 0: the weak coupling expansion 
of the pressure is organized in single powers (and logs) of g. The relevant effective theory is given 
by dimensionally reduced electrostatic QCD (EQCD) which for the pressure has been worked out 
up to but not including order g^ which is where nonperturbative physics from magnetostatic QCD 
(MQCD) starts to contribute. This barrier is indicated by a thick red line in Fig. 1171 

Because g is treated as an arbitrarily small parameter, everything in Fig. El with the exception 
of the border at x = corresponds to the regime T <C /i, namely T ~ g^fi with x > 0. As long as 
X < 1, T is parametrically larger than the Debye mass ~ gn, and so dimensional reduction should 
still be applicable. However, each coefficient of the original series at x = now has to be expanded 
in powers of T/fi ~ g^ . The 2-loop pressure contribution for example yields three different terms 
for X > 0: one is proportional to fi'^ and thus is always of order g'^, another — proportional to fi^T"^ 
— gives the line y = 2 -|- 2x and the third term proportional to T^ produces the line y = 2 -\- Ax. 
Starting with the plasmon term which is of order ^^ at x = 0, we obtain an infinite series of 
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FIG. 17: The structure of the weak-couphng expansion of the interaction pressure p — psB at parametricaUy 
small r//i as a function of the power a; in T '--^ 3^M- At T ^ fi, i.e. x — 0, the expansion involves orders 2, 3, 
4,. . . in g (logarithms in g are not made explicit); at T ~ 5^' ^M where dimensional reduction overlaps with 
HTL/HDL resummation, the series in g involves powers 2, 3, |, 4, |,- • ■ ! at T '--^ g^i, where dimensional 
reduction ceases to be applicable, the expansion is again in even powers of g (and logs) with coefficients that 
at even smaller temperatures can be expanded in a series involving fractional powers of T (beginning with 
2, |, ^, 4, ^,. . . ) and corresponding powers 2 + 2x, 2 + |x, ... of g. While subleading in the pressure, the 
latter contributions give the leading-order "anomalous" (non-Fermi-liquid) contributions to the interaction 
part of the entropy and specific heat at low temperature. Existing results for the various contributions are 
represented by full lines, as-yet undetermined contributions by dashed and dash-dotted lines. The non- 
perturbative barrier from the scale of magnetostatic confinement (magnetic screening mass) is indicated 
by the thick red line marked "non-perturbative" . The region below it and up to a; = 1 is the regime of 
electrostatic QCD (EQCD), while for a; > 1 the relevant effective theory is given by non-static hard dense 
loops (HDL). 



higher-order terms for x > 0. These arise from the expansion of the third power of the Debye 
mass parameter in powers of T/fi, and, for subsequent terms in the dimensional reduction result, 
also from the expansion of the special functions H(n, z). Because both the Debye mass and the b^ 
functions can be expanded in even powers of T/fi, the lines emanating from their starting points 
at X = come with slopes differing by two units. The terms proportional to g^ and ^r^ at x = 
involve a single overall power of T, so the lines emanating from these have slopes 1, 3, 5, . . . , 
whereas the term proportional to g'^ (or gHng) has also T-independent parts and thus gives rise 
to lines with slopes 0, 2, 4, ... . In Eq. (|2.6j) we have seen how this gives rise to a new series in g 
at a; = 1/2, and Fig. [T7I illustrates how the individual terms of order 2, 3, |, 4, |, . . . are produced 
from the various coefficients of the expansion at x = 0. 

Moving on to the border of applicability of the dimensional reduction results, x = 1, we see that 
all lines converge to points corresponding to an expansion in even powers of g (and also involving 
Ing). As noted before, for x > 1 the relevant effective theory is the one given by non-static hard 
dense loops. Their resummation is necessary to obtain the classic Freedman-McLerran (FMcL) 
result to order g'^ (again accompanied by a logarithmic term) at T = as well as the leading 
thermal corrections to the interaction pressure. In a low-T expansion these T-dependent terms 
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start with a contribution of order g'^T'^fi^ln{T/gfi) and then involve fractional powers T^'^, T^^'^, 
T^lnT, T^'^'^, . . . such that the corresponding lines in Fig. El (labeled by the exponent of T) 
meet at x = 1 and effective order g'^. At this point, the leading T-dependent contributions are of 
the same order as the three-loop T = (FMcL) pressure contribution and remain more important 
than the undetermined four-loop T = term even for parametrically lower temperatures as long 
as X < 2 (i.e. T >> g'^fi). For the entropy and specific heat, for which the zero-temperature 
contribution to the pressure drops out, these T-dependent terms represent the leading interaction 
contributions down to arbitrarily low temperatures. The TlnT behavior of the entropy (as well 
as of the specific heat) is characterized by "anomalous" non-Fermi-liquid behavior, caused by the 
only weakly (dynamically) screened quasi-static magnetic interactions with an effective frequency- 
dependent screening mass, displayed in Eq. (|2.9j) . 

As suggested by Fig. Eland shown in detail in the previous section, the HDL-resummed thermal 
pressure contributions responsible for the non-Fermi-liquid behavior at T <^ gfj, match smoothly 
to the perturbative effects at T ^ gfi described by EQCD. As the temperature is increased, elec- 
trostatic screening replaces dynamical magnetic screening as the dominant collective phenomenon 
also in the T-dependent contributions. For T parametrically larger than gfj, (i.e. x < 1) the resum- 
mation of HDL self energies needs to be trivially extended to HTL self energies to avoid accuracy 
loss. When added to the zero temperature 0{g^) result, this gives an expression that gives the 
pressure for all temperatures and chemical potentials up to an error of order gr™™(''+2^'6) (or g'^^^ 
throughout in the case of the entropy, for which the unknown four-loop T=0 pressure drops out). 

From the "flow" of the various perturbative contributions as a function of x in Fig. 1171 one 
notices that a single expression aiming to be valid both for x > 1 and x < 1 needs to keep track 
of contributions which are perhaps higher-order and irrelevant in some region but essential in 
another. The novel approach we have presented here does so by resumming the complete one- 
loop gluon self-energy in all IR sensitive graphs while treating the infrared-safe 2GI diagrams 
perturbatively. To the extent that we have worked it out, this procedure covers both x > 1 and 
X < 1 with an error of order (^min(5+x,6) -^j^ici^ improves over the HTL/HDL result in the region 
X < 1 by including the contributions of all relevant three- loop graphs. A drawback compared to the 
HTL/HDL resummation schemes is however that the resummation of the complete gluon self-energy 
leads to gauge-dependent higher-order contributions whose unphysical nature is highlighted by the 
appearance of spacelike poles in the logarithmic resummation integrand with momenta ~ g^T and 
also of an unphysical damping constant (with an incorrect sign) oc g^T. For our expression for the 
pressure, the effect of these problems is, however, only of the order of the nonperturbative MQCD 
contributions, i.e. g^, so it has not hindered us from confirming and thus validating the results 
obtained through dimensional reduction or the HTL/HDL approach. 

VI. CONCLUSIONS AND OUTLOOK 

In this paper, we have constructed a novel resummation scheme designed to reproduce the weak 
coupling expansion of the QCD pressure up to order g^ on the entire fi-T plane. We have used it to 
provide an independent check of practically all existing perturbative results. In particular, we have 
performed the first explicit test on the validity of dimensional reduction for values of fj,/T far beyond 
the capability of present-day lattice techniques, thus verifying that dimensionally reduced effective 
theories provide a solid description of the perturbative physics up to in principle arbitrarily large 
values of fi/T as long as vrT > ttiq. At temperatures parametrically smaller than the chemical 
potential, we have on the other hand reproduced numerically all the results of the HTL/HDL 
resummation schemes, verifying their validity and highlighting the smooth transition taking place 
in the perturbative expansion of the pressure as one moves from the region of dimensional reduction 



33 



towards the zero-temperature limit. 

Based on our numerical results from Section HVl the dimensional reduction result for the QCD 
pressure appears to be provide a remarkably good approximation for this quantity down to the 
point where the T-dependent contribution to the interaction pressure, Sp, ceases to be negative 
(cf. Figs. [Tp) which happens at T ~ 0.2mD. Since the dimensional reduction result to order g^lng 
combined with optimized choices of the renormalization scale has turned out to agree rather well 
with lattice results, both at zero chemical potential [3, |4l| and for /i ~ T" [3, |42|, our present 
findings in fact suggest a remarkably wide practical range of applicability for the dimensional 
reduction method and results. 

Progressing down on the temperature axis to T < 0.2777-d, one eventually has to switch to the 
nonstatic resummation schemes provided either by our new approach or by the calculationally 
much simpler HTL resummation of Eq. (|2.15|) . At such low temperatures, the pressure can — up 
to but not including order g^ — be approximated by the Preedman-McLerran result plus positive 
contributions from the Stefan-Boltzmann terms as well as the interaction pressure 5p. The latter 
of these is the source of the non-Fermi-liquid behavior of the entropy and specific heat. 

While we believe to have thoroughly clarified the nature of perturbative expansions of the pres- 
sure in different regimes of the fi-T plane, our new approach is, as of today, yet to produce results 
for the pressure beyond what has already been achieved through either dimensional reduction at 
X < 1, the HTL/HDL resummation schemes at x > 1 or the Freedman-McLerran result at T = 0. 
Its present relative error of order gr™™(5+2^,6) ^^j^ jj^ principle be reduced through the inclusion of 
the two-loop gluon polarization tensor into the resummation of the ring diagrams and in addition 
by taking the contributions of non-static modes into account in the multiple sums of Fig. 2.b- 
d. For the pressure, this would bring the accuracy of our new approach up to the one currently 
achieved by dimensional reduction calculations (excluding the already known 0{g^lngT'^{T'^ + fi'^)) 
term), so that the error, up to logarithms, would be uniformly (for all values of x) of order g^lng, 
corresponding to the line marked "4-loop T = pressure" in Fig. 1171 This would then unify all 
existing perturbative results for the pressure of QCD, while for the entropy it would moreover 
lead to genuinely new results. Apart from increasing e.g. the accuracy of the entropy result at 
X = 1/2 to order g^^'"^ (green open dots in Fig. 117(1 ''^^. it would push the error in the T-dependent 
part of the pressure up to the line denoted in Fig. El by "4-loop T contribution" and thus, for 
X > 1, include the so far unknown order g^^i^T corrections to the non- Fermi-liquid terms in the 
entropy and the specific heat. Considering the difficulties caused by the gauge-dependent parts of 
the gluon self-energy, it seems that such an extension should probably aim at keeping only gauge- 
independent contributions such as HTL self energies in the ring diagrams and treating corrections 
to those self-energies in a perturbative manner. Work towards this goal is currently in progress. 
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^^ The highest purely perturbatively calculable order at x = 1/2 is g^^'^ which would require a calculation of the 
contributions of order g^fi^T^ and g^ fJ.'^T for the pressure. 
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APPENDIX A: DIMENSIONAL REDUCTION RESULT AT FINITE T AND n 

To order g^ and simplified by assuming equal chemical potentials, the result of dimensional 
reduction for the QCD pressure at finite T and //, obtained in Ref. [a], reads 
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APPENDIX B: THE ANALYTIC VALUES OF VARIOUS PARTS OF THE PRESSURE 

In this Appendix, we collect some calculational details on the different pieces of the analytical 
part of the pressure, Pani, as defined in ()3.17|1 . The sum of the 2GI diagrams was already given in 
Eq. (|3.1j) and the piece Pb in Eq. (|3.16|) . 

1. The Vac- Vac and Vac-Mat diagrams 

The evaluation of the two special diagrams of Fig. 2, dubbed Vac- Vac and Vac-Mat based on 
their self energy insertions, is relatively straightforward. Inserting the form of Ii'f^^{P)\^^c from 
Eq. (|C4|) into the Feynman gauge expressions of the graphs, contracting all Lorentz and color 
indices and taking into account the symmetry factors 1/4 and 1/2, respectively, we readily obtain 
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where we in the latter case have two two-loop sum-integrals to evaluate. These can be taken care 
of using methods developed in Refs. [J, 1^, which leads to the final result 
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2. The ordinary ring sum 

The UV subtraction term p^^^ of the ordinary ring sum, defined in Eq. 1)3. 8j) . is straightforward 
to evaluate and produces the result 
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With the IR regulating term p]^^^, we find it most convenient to first transform the Matsubara sum 
into a contour integral in the standard way. After the term proportional to a bosonic distribution 
function is analytically continued to Minkowski space (see Appendix D), this leads to the result 
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The sum of the first two terms is evaluated with standard methods and gives 
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while the third term (where one can set e = 0) is left in its present form until Appendix D where it 
is seen to explicitly cancel with certain parts of p^^^". Modifications to plf^ induced by a magnetic 
mass are discussed after Eq. ()D5|1 . 



3. Small and large fi/T limits of pani 

For convenience, we provide here formulae for various limits of the function p^^i of Eq. H3.2U() . 
derived using the results of Appendix D of Ref. ^S] for the behavior of the H functions. First, in 
the limit // ^ we obtain through an expansion in fl = ^^ 
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which contains e.g. the CpTp dependent part of the hnear quark number susceptibihties found in 
Ref. ^44j. Taking the opposite hmit, T — > 0, we on the other hand get 
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This yields an entirely finite expression at T = 0, where the scale dependence of the result naturally 
coincides with that of Eq. (|2.7|) . 

APPENDIX C: PROPERTIES OF THE ONE-LOOP GLUON POLARIZATION TENSOR 



For reference, we note that the one-loop gluon polarization tensor in the Feynman gauge can 
be written in the form 
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The vacuum (T, /i ^ 0) limit of the above formula is uniquely determined by replacing the sum- 
integrals by 4 — 2e-dimensional integrals which produces 
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As all our analytic formulae, the above expressions have been given in Euclidean space. For the 
numerical evaluation of pl'^^ it is, however, convenient to also reproduce the known formulae for 
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the transverse and longitudinal pieces of the vacuum-subtracted self energy in Minkowski space. 
Following Ref. |4f)l |. we obtain 

nL((?o,9) = -g^% {2TFHf + CaH^) , (C6) 
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where the different functions read 
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These expressions are valid for all complex go in a rotation from Euclidean space qo = iuj to 
Minkowski space q^ = lo + ie with e > 0. For the analytic continuation into the region with e < 0, 
see Ref. 0. 

The corresponding HTL and HDL expressions for the above functions can be extracted by 
demanding that the external momenta go and q be small compared to the temperature T and/or 
the chemical potential /i. This amounts to expanding the integrand in inverse powers of k to 
leading order, since the main contribution is expected to come from large loop momenta k. The 
relevant physical scale is given by the Debye mass m^. In this limit, the equations HC8|) to HC11|1 
reduce to 

/ff- = ?|l + |llogf^i±^H. (C17) 
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The HDL expressions are obtained by setting T = here which makes both G™^ and Hf^^ vanish. 

APPENDIX D: NUMERICAL EVALUATION OF p^^^'' 

The numerical evaluation of the sum-integrals introduced in Sec. HI is performed by converting 
the sums over Matsubara frequencies into contour integrals where the integrand is multiplied by 
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cot(a;„/2T) and the contour encircles the poles of this function (which lie on the real axis at 
cOn = 2TTnT and have the residue 2T). We then change the integration variable from qq = iujn, 
with real ujn, to qo + ie, with real qo and e > 0, deforming the integration contour in the usual way 
(for more details on this, see e.g. Ref. [^d])- Through this procedure, we obtain the generic result 
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^dqo 
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where we have used the identity coth ^ = 1 + 2nh. This is valid for any function / satisfying 
Hqo + ie) = f*iQo — ie) = fi—Qo — ie) = /*(— 90 + ie), provided that the great arc contribution 
to the integral vanishes and there are no other poles or branch cuts between the Euclidean and 
Minkowskian axes. 

After these standard manipulations, we split the above integral into the UV safe "rif, part" 
which we calculate in Minkowski space, i.e. along the real qq axis, and the "non-n^ part" that we 
rotate back to Euclidean space^^ before integration through 

"^Im/(go + ^6) - T— Re/(H- (D2) 

^ Jo ^ 

In this way, we take full advantage of the UV cutoff properties of rif, in Minkowski space while 
the possibly UV problematic non-n^ parts are treated respecting the Euclidean symmetry. Further 
discussion on the reasons for this separation can be found from Refs. |l5l.l32]. 

Specializing to the evaluation of the function pfi^g") there is a further subtlety in the non-n^ 
part that we have to take into account. As discussed in Sec. 1111 Fl there is an unphysical pole 
in the integrand of Eq. H3.6() which implies that we cannot deform our integration path from the 
Euclidean to the Minkowskian axis quite as suggested by Eq. HD1|) . Instead, we have to try to avoid 
numerically dangerous singularities on the complex qo plane by encircling them at a safe distance. 
This aspect of our calculation is discussed in more detail below in ADDendix lD2al 

After all the separations, we can write the function p^-!^^" as the sum of several individual 
contributions according to 

p^::: = pf + pf + < + p^B + pt + ptd + pZ- (d3) 

Here, pf and pf denote two distinct parts of the Euclidean integral that are separated by a mo- 
mentum cut-off A. The pieces p^ and p^ on the other hand stand for longitudinal and transverse 
Minkowskian contributions to p^-^^^" integrated along specific segments of the real qo-axis and p^D 
and p^Q for the corresponding quantities integrated along deformed paths (see Fig. ()D2 a|) ). In 
addition, m refers to corrections due to the IR regulator mass introduced in Sec. Ill C that can be 
explicitly extracted from the rest. Each of these terms will be discussed in detail in the following. 

1. Euclidean space contributions 

For practical numerical reasons, we divide the Euclidean space "non-nfe" contribution to Eq. ()D3|1 
into two parts pf + pf by introducing a four-momentum cut-off A and having the first piece 



^® The quantities evaluated in Minkowskian metric (+,—,—,—) (implying Q^ ~ qo — q^) are denoted by an upper 
index M, while Euclidean pieces, where qo is analytically continued back to Euclidean space (go —> —iqo = <^), are 
indicated by the use of the "Matsubara variable" u and the upper index E. We have dropped the index n from cOn 
as it now has become a continuous variable. 
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correspond to the contribution of momenta with \Q\ < A. For this, we obtain by first performing 
the integration over an Euchdean four-sphere 
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with u! = —iqo = Q cos a and q = Q sin a. The integrand pf^^iuj, q) is given by 
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where the magnetic mass m^^^, introduced in Eq. (|l-{.24|) . is needed in the transverse part in order 
to prevent a negative argument of the logarithm as w^ + g^ -^ 0. This provides a cutoff for IR 
divergences at the scale m^^g = cj-^g'^CAT, with c/ > 1 being a "magnetic factor" that we can 
vary to verify that the dependence of the result on it is beyond 0{g'^). 

In addition to regulating IR divergencies, the inclusion of the magnetic mass alters the UV 
behavior of the first term of the above integrand which is why we included it also in the coefficient of 
the UV regulating last term of the first line of Eq. ()D5|) . In order to keep the final result independent 
of the arbitrary mass parameter m, we have to take the magnetic mass m,-,-,^^ into account also in 
the calculation of the function p™^^ of Eqs. (|3.7|) and ()B5|) . which amounts to replacing C^ in those 
equations as well as later in Eqs. ()D10|) and (|D15|) by C'^ — 8g^CArn^^^^^ + 24:g^rn^^^^, as shown below 
in Eq. (lD7|) . 



a. High momentum expansion in Euclidean space 



To obtain the high-momentum Euclidean contribution p^ to Eq. HD3|) in the most effective way, 
we first integrate over a, then expand 11^ and IIl in the limit of large momenta (which is possible 
analytically order by order in Q/T) and finally perform the integration over Q. This produces a 
series of the general form 
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Phigh.O ~r ^high,2 ' Phigh,4 I ••• 



(D6) 



where the lowest orders are given by 
d-A 



-Phigh,0 



Phigh,2 



7687r2 
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g^T* [C\ - Sg^CArnl^^ + 24/mi,J log 



34567r2A2 



g^T' 



' niag 



9^C^ 



2^2 



A2 



m 



m- 



+ A2 



m- 



+ A2 



^3 



.6— -6 



36g''CArK.., + 72g'^K.., 



^ + 8Ca + 14T^ ^ ^gj.^ _2 ^ gg^^ _4 



(D7) 



(D8) 



with m^iag = "imag/(5^r) = Cf-^CA (with Cf > 1) and // = ^/{2-kT) as before. In our numerical 
implementation, we included terms up to Phigh,io and used a cut-off A ~ 160 y^^ + /?, for which the 
error in the sum of pf and p^ due to the introduction of the cut-off and the use of a high-momentum 
expansion in the second part is negligible. 
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FIG. 18: Symbolic illustration of the analytic structure of the integrand oi p^i^g" on the complex plane. The 
unphysical pole is avoided by using the complex integration path described in the text. 

b. IR mass in Euclidean space 

The IR regulating mass m only appears in Eq. ()3.6|) in two terms that are independent of the 
self energies, and thus the difference of p^"'*" evaluated with masses mo and m is given by 



p^::^^-p^::T(^o) 



2 „4rr4 
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p[{P^+m^f (p2+^2)2 



(D9) 



From Eq. HiJ6() . we see that the Euclidean contribution to this sum- integral is 



2 „4Tn4 



p^:;;fM-p^::f("^o) 



dACigy- 

24(47r) 



log- 



m 



mo 



(DIO) 



which implies that the cancelation of m in the sum of the Euclidean contributions to p^f^" can be 
tested independently of the Minkowskian ones by verifying that the m dependence of the numerical 
result for pf +P2 is exactly of the above form. We have done so with the expected positive result. 

2. Minkowski space contributions 

a. Contour deformation 

As discussed in Sec. 111111 even after the introduction of the magnetic mass in p^.^^" there is 
an unphysical pole remaining on the complex qq plane which prohibits the use of the preferred 
"nfe" integration path along the Minkowskian axis. Instead, we have to choose our integration con- 
tour so that we avoid the unphysical pole, illustrated in Fig. ^J Since a standard path along the 
Minkowskian axis would fail at the branch cut from qq = q to the unphysical pole, our numerical 
integration contour goes around the problematic zones at a safe distance, so that (given the van- 
ishing great arc contribution) the final integration is equivalent to integrating along the Euclidean 



42 





1 












^^^"""^.To 


r^ 








/ 


0.5 










/ 


0.25 


.^►j- 








1-1 


-0.5 

-0.25 

-0.5 

'~'"^^^^^) . 7 5 




0.5 


1 


1.5 2 




-1 







FIG. 19: The shape of the Debye mass cutoff and the numerical integration contour chosen on the complex 
go plane, in units of r. The full line shows the deformed integration path along complex go while the dotted 
line indicates usual Minkowski space integration, extending to go ^ oo. The dash-dotted line shows the 
integration path if it is restricted to stay in Minkowski space up to go = q/S. The parameters in this and 
the following two Figures are chosen (a) for the big circle shape (g ~ .5r) 



axis — or just summing the discrete Matsubara frequencies we started from. The validity of this 
procedure can be checked afterwards by varying the safety distance between our integration path 
and the branch cuts. 

The unphysical poles are expected to appear in the vicinity of the lightcone. We thus want our 
integration path to avoid the region \q^ — (^\ <r^ where r is chosen to be of the order of the Debye 
mass r = fcrriu and /c > 1 is an arbitrary cutoff factor on which the final result should not depend. 
We will, however, not integrate exactly along this path, but instead choose a rectangular encasing 
boundary. To determine the shape of this optimal rectangular path, we write the condition 



kl 



\ia + ibf 



q \ < r^ 



in terms of the real quantities a = Re qq and b = Im qq which gives the relation 



Aa'^b'^ + 



i'? < r\ 



(Dll) 



(D12) 



Solving for the extremal points for a and b from here and denoting x = ^Jr"^ + q^ , 
z = r^/(2g), we find the three cases 



= X^ 



V2 



r 
< a < r 



> r 



qo G {0, iy, x + iy, x} (big circle) 

go S {0, iz^ X + iz, x} (eyeglasses) 
9o G {Vi y + iz,x + iz, x} (two circles) 



where the notation qq £ {u, v, w, z} implies that the integration path on the complex qq plane pro- 
ceeds along straight lines through these points. In practise, it is convenient to start the integration 
along the Minkowskian axis also in the two first cases, so from to q/2 the integration is always 
performed in Minkowski space. The three contours are depicted in Figs. ITOl- 1211 

With the above integration paths, the notation in Eq. (|D3|1 becomes obvious. The pieces p^ 
and p^ denote the integrations performed along the Minkowskian axis (dashed lines), making use 
of optimizations [S^], while p^o ^^id p^^ are computed using the deformed paths (full lines or 



dash-dotted lines) for complex go- The Minkowskian contributions p^ and p^ are calculated as 



Ph/T 



'27r3 



dqolniipttiQo,q)) 



(D13) 
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FIG. 20: (b) for the eyeglasses shape (g = .97r) 
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FIG. 21: and (c) for the two circles shape (g = 1.2r) 



where go,min denotes the starting point of the green dashed lines in Figs. El~l^ (with the exception 
of the two circles shape for which the integration contour is composed of two disjoint pieces) and 
p?^j is the straightforward analytic continuation of Eq. I)D5|1 from Euclidean to Minkowski space. In 
Minkowski space, the Boltzmann factor n^ provides a natural UV cutoff and the IR mass dependent 
term can be evaluated separately (see below). 



The contribution of the deformed paths, pjfn +Ptd' i^ computed similarly as 



/^LD/TD 



dA_ 

2^3 



TTZs / ^^^^ 



dqolm{p^^^{qo,q)). 



(D14) 
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b. IR mass in Minkowski space 



The so far unevaluated last contribution to Eq. (|D3|) is p^ which is the IR mass correction term 
corresponding to the m-dependent terms of Eq. (|3.6() . This integral can be computed separately 
from the other Minkowskian pieces, since the bosonic distribution function provides a natural UV 
cutoff for the Minkowski space calculations, making the subtraction of the quartic term from pf;"'*" 
unnecessary. Writing p^ down explicitly (see the discussion below Eq. (|D5|1 on how to include 



m„ 



we obtain 



2 ^4rr4 



Pn 



24 



dqo f d?q 



TT 
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n^Im 



{Q' 



m" 



(D15) 



which, not surprisingly, is observed to exactly cancel the third term of Eq. HBSf) . Thus we see that 
all dependence on the regulator mass m cancels between the different Minkowskian contributions 
to plll° which together with our observation after Eq. (|D10|) explicitly verifies the independence of 
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our final result piy of m. 



[1] E. V. Shuryak, Sov. Phys. JETP 47 (1978) 212 [Zli. Eksp. Teor. Fiz. 74 (1978) 408]. S. A. Chin, Pliys. 

Lett. B 78 (1978) 552. 
[2] J. I. Kapusta, Nucl. Pliys. B 148 (1979) 461. 
[3] T. Toimela, Phys. Lett. B 124 (1983) 407. 



[4] P. Arnold and C. X. Zhai, Phys. Rev. D 50 (1994) 7603 |hep-ph/940 82m ; Phys. Rev. D 51 (1995) 
1906 hep-ph/9410360 . 

[5] C. X. Zhai and B. Kastening, Phys. Rev. D 52 (1995) 7 232 |hep-ph/9507380| . 

[6] E. Braaten and A. Nieto, Phys. Rev. D 53 (1996) 3421 hep-ph/9510408 . 

[7] K. Kajantie, M . Laine, K. Rummukainen and Y. Schroder, Phys. Rev. D 67 (2003) 105008 
|hep-ph/0211321| . 



[8] A. Vuorinen, Phys. Rev. D 68 (2003) 054017 hep-ph/0305183 . 

[9] B. A. Freedman and L. D. McLerran, Phys. Rev. D 16 (1977) 1169; V. Baluni, Phys. Rev. D 17 (1978) 

2092. 

[10] A. Ipp, A. Gerhold and A. Rebhan, Phys. Rev. D 69 (2004) 011901 |hep-ph/03090T9 . 

[11] A. Roberge and N. Weiss, Nucl. Phys. B 275 (1986) 734. 

[12] J. C. Taylor and S. M. H. Wong, Nucl. Phys. B 346, 115 (1990); E. Braaten and R. D. Pisarski, Phys. 

Rev. D 45, 1827 (1992). 

[13] D. T. Son, Phys. Rev. D 59 (1999) 094019 |hep-ph/981 2287|. 

[14] A. Gynther, Phys. Rev. D 68 (2003) 016001 |hep-^V 03036l9]. 

[15] A. Ipp and A. Rebhan, JHEP 0306 (2003) 032 ■hep-ph/0305030'. 

[16] A. Gerhold, A. Ipp and A. Rebhan, Phys. Rev. D 70 (2004) 105015 hep-ph/040608 7 . 

[17] A. Vuorinen, "The pressure of QCD at finite temperature and quark number density," PhD thesis 

hep-ph/0402242". 
[18] M. Laine and Y. Schroder, "Quark mass thresholds in QCD ther modynamics," |hep-ph/0603048, , 
[19] E. S. Fraga and P. Romatschke, Phys. Rev. D 71, 105014 (2005) |hep-ph/0 412298 . ' 
[20] A. Hart, M. Laine and O. Philipsen, Nucl. Phys. B 586, 443 (2000) ■hep-ph/0004060'. 



[21] C. P. Korthals Ahes, R. D. Pisarski and A. Sinkovics, Phys. Rev. D 61, 056007 (20 00) hep-ph/9904305| . 



[22] J. P. Blaizot, E. lancu and A. Rebhan, Phys. Lett. B 523, 143 (2001) hep-ph/ 0110369 

[23] T. Holstein, R. E. Norton, and P. Pincus, Phys. Rev. B8, 2649 (1973). 

[24] A. Gerhold and A. Rebhan, Phys. Rev. D 71, 085010 (2005) hep-p h/0501089 . 

[25] J . O. Andersen, E. Braaten, E. Petitgirard, and M. Strickland, Phys. Rev. D66, 085016 (2002) 

|hep-ph/0205085' ; J. O. Andersen, E. Petitgirard, and M. Strickland, Phys. Rev. D70, 045001 (2004) 

'hep-ph/0302069 . 
[26] J. P. Blaizot, E. lancu, and A. Rebhan, Phys. Rev. Lett. 83, 2906 (1999)_jhep-ph/9906340 ; Phys. Lett. 

B470, 181 (1999) hep-ph/9910309 ; Phys. Rev. D63, 065003 (2001) 'hep^ph/0005003 . 
[27] E. S. Fraga, R. D. Pisarski and J. Schaffner-Bielich, Phys. Rev. D 63, 121702 (2001) hep-ph/0101143| . 
[28] A. D. Linde, Phys. Lett. B 96, 289 (1980). 
[29] O. K. Kalashnikov and V. V. Klimov, Sov. J. Nucl. Phys. 33, 443 (1981) [Yad. Fiz. 33, 848 (1980)]. 

[30] O. K. Kalashnikov, Fortsch. Phys. 32, 525 (1984). 

[31] U. Kraemmer and A. Rebhan, Rept. Prog. Phys. 67, 351 (2004) |hep-ph/0310337 . 



[32] G. D. Moore, JHEP 0210 (2002) 055 |hep-ph/0209190| ; A. Ipp, G. D. Moore and A. Rebhan, JHEP 

0301 (2003) 037 hep-ph/0301057 . 

[33] J. Kapusta, "Finite temperature field theory", Cambridge University Press (1989). 

[34] O. K. Kalashnikov and V. V. Klimov, Sov. J. Nucl. Phys. 31, 699 (1980) [Yad. Fiz. 31, 1357 (1980)]. 

[35] K. Kajantie and J. I. Kapusta, Annals Phys. 160, 477 (1985). 

[36] U. W. Heinz, K. Kajantie and T. Toimela, Annals Phys. 176, 218 (1987). 

[37] E. Braaten and R. D. Pisarski, Phys. Rev. D 42, 2156 (1990). 

[38] T. Schafer and K. Schwenzer, Phys. Rev. D 70, 0540 07 (2004) !ar Xiv:hep-ph/0405053 . 

[39] A. Schmitt, Phys. Rev. D 71, 054016 (2005) |nucl-th/0412033 , 



[40] W. E. Brown, J. T. Liu and H. c. Ren, Phys. Rev. D 61, 114012 (2000) |arXiv:hep-ph/9908248| ; 



45 



Phys. Rev. D 62, 054013, 054016 (2000) arXiv:hep-ph/9912409,0003199]; A. Schmitt, Q. Wang and 
D. H. Rischkc, Phys. Rev. D 66, 114010 (2002) arXiv:nucl-th/0209050 . 

[41] J. P. Blaizot, E. lancu and A. Rebhan, Phys. Rev. D 68, 025011 (2003) he p-ph/0303045 |. 

[42] A. Ipp, A. Rebhan and A. Vuorinen, Phys. Rev. D 69, 077901 (2004) hep- ph/0311200| . 

[43] J. P. Blaizot, A. Ipp and A. Rebhan, "Study of the gluon propagator in the \arge-Nf Hmit at fi- 
nite temperature and chemical potential for weak and strong couplings," Annals of Physics (in print) 
|hep-ph/0508317 . 

[44] A. Vuorinen, Phys. Rev. D 67 (2003) 074032 |hep-ph/0212283 . 

[45] H. A. Weldon, Phys. Rev. D 26, 1394 (1982). 



